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CHAPTER  I 
INTRODUCTION 

The  frequency  of  discovery  of  gas  condensate  reservoirs, 
as  well  as  their  importance  to  the  industry,  has  increased 
greatly  in  recent  years  with  the  trend  toward  drilling  into 
the  deeper  environment  of  higher  temperatures  and  pressures. 
Craft  and  Hawkins   point  out  well  test  data  confirm  this 
trend  as  shown  in  Figure  1.   This  graph  shows  the  discovery 
trend  for  17  parishes  in  southwest  Louisiana.   The  reservoirs 
were  separated  into  oil  and  gas  or  ga s - co nde ns a t e  types  on 
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Fig.   1.    Discovery  frequency  of 
oil  and  gas  or  ga s - c onde ns a t e 
reservoirs  versus  depth.    (Craft 
and  Hawkins,  Applied  Petroleum 
Eng  i  n eer i  n g .  ) 


the  basis  of  well  test  gas  oil  ratios  and  the  API  gravity  of 
the  produced  liquid.   While  oil  discoveries  predominate  at 
depths  less  than  8000  feet,  gas  and  gas  - condens a te  discover- 
ies predominate  below  10,000  feet.   The  decline  in  discover- 
ies below  12,000  feet  is  due  to  the  fewer  number  of  wells 
drilled  below  that  depth  rather  than  a  drop  in  the  occurrence 
of  hydrocarbons. 

Many  of  these  high  pressure  condensate  reservoirs  are 
discovered  in  conditions  above  their  dew  points.   With  pres- 
sure reduction,  when  placed  on  production,  they  allow  retro- 
grade condensation  of  the  valuable  heavier  hydrocarbons  within 
the  reservoir  rock.   Once  condensation  has  taken  place  within 
the  reservoir,  this  liquid  cannot  be  produced  without  con- 
siderable expense.   The  mechanisms  by  which  this  condensation 
may  be  retrieved  is  either  by  r evapor i za t ion  or  by  repressuri- 
zation  by  injected  dry  gas.   However,  several  pore  volumes  of 
dry  gas  must  usually  be  injected  before  sufficient  revaporiza- 
tion  occurs,  making  this  process  costly.   Even  then,  the  dry 
gas  may  not  sweep  all  the  volume  in  which  condensation  has 
taken  place.   In  order  to  retrieve  this  condensate,  which  is 
many  times  more  valuable  than  the  gas  itself,  or  to  prevent 
its  condensation  within  the  reservoir,  it  is  necessary  to  in- 
vestigate the  possibility  of  dry  gas  cycling  and  pressure 
ma  intenance. 

Cycling  of  the  dry  residue  gas  may  act  to  maintain  the 
pressure  in  the  reservoir  at  or  near  the  dew  point,  thereby 


preventing  the  occurrence  of  the  retrograde  condensation 
phenomena.   It  may  also  tend  to  revaporize  the  condensate  and 
sweep  it,  along  with  the  cycled  gas,  out  of  the  reservoir. 

Even  though  it  is  apparent  that  cycling  a  reservoir 
which  is  susceptible  to  the  retrograde  phenomena  will  produce 
more  of  the  valuable  hydrocarbons  in  place,  the  economic  fac- 
tors involved  are  all  important  in  deciding  whether  to  invest 
in  this  procedure.   Roughly,  a  reservoir  that  contains  as  much 
as  fifty  barrels  of  condensate  per  million  standard  cubic  feet 
of  gas  is  a  good  prospect  for  gas  cycling.   On  the  one  hand, 
the  dry  gas  sales  will  be  postponed  until  the  cycling  is  prac- 
tically complete,  resulting  in  decreased  present  value  of  the 
gas.   However,  it  is  generally  possible  to  produce  the  con- 
densate faster,  hence  earlier,  when  cycling  is  adopted  and 
more  dollars  are  generally  obtained  due  to  the  increased  ulti- 
mate recovery.   The  purpose  of  this  thesis  is  to  evaluate  whether 
or  not  these  advantages  are  worth  the  greater  investment  in- 
volved for  the  field  studied. 

Data  are  obtained  from  an  actual  field  on  the  Gulf  Coast 
of  Louisiana  containing  approximately  600  billion  SCF  of  gas. 
The  size  of  the  reservoir  was  determined  from  geologic  data 
and  from  production  data  to  date.   Future  prediction  of  reser- 
voir performance  was  made  for  pressure  depletion,  normal  cycling 
and  pressure  maintenance  methods  of  exploitation.   An  economic 
comparison  of  the  cases  studied  is  presented. 


CHAPTER  II 
DISCUSSION 


A .  Nature  of  Gas  Condensate  Reservoirs 

Gas  reservoirs  are  defined  as  having  a  lean  gas,  con- 
taining a  minimum  amount  of  the  heavier  hydrocarbons,  where 
the  composition  of  this  produced  gas  does  not  change  with 
pressure  reduction  as  the  reservoir  is  depleted.   Oil  reser- 
voirs range  in  gas  content  from  zero,  or  dead  oil,  to  a  gas- 
oil  ratio  of  a  few  thousand  cubic  feet  per  barrel.   In  gas 
reservoirs,  on  the  other  hand,  the  liquid  may  be  vaporized 
into  100,000  cu.  ft.  per  barrel  or  more.   Gas  condensate 
reservoirs  are  classified  between  these  two  extremes.   Gas 
condensate  reservoirs,  or  distillate  reservoirs  as  they  are 
sometimes  called,  usually  produce  light  colored  or  colorless 
stock  tank  liquids  with  gravities  above  45  degrees  A.P.I,  at 
gas-oil  ratios  in  the  range  of  3,000-100,000  SCF/bbl.   The 
production  from  these  reservoirs  is  predominately  gas  from 
which  some  liquid  is  condensed  in  the  surface  separators, 
which  phenomenon  gives  them  the  name  gas  - co ndens a t e . 

Another  definition  of  the  three  types  of  fields  described 
can  be  derived  from  an  interpretation  of  the  pr es sur e- t emp era tur e 
phase  diagrams.    (See  Figure  2  for  a  typical  diagram.)   When- 
ever a  reservoir  fluid  is  discovered  at  a  temperature  above  the 


c r ic onden therm,  or  maximum  two  phase  temperature  for  the  fluid, 
the  composition  of  the  produced  fluid  will  not  change  as  the 
reservoir  is  depleted  since  the  reservoir  temperature  can  be 
expected  to  remain  fairly  constant  (Path  A-A   on  the  diagram). 
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Fig.  2.   Pr es sur e- temp  era tur e  phase 
diagram  of  a  reservoir  fluid.    (Craft 
and  Hawkins,  Applied  Petroleum  Rese r  v  o  i  r 
Engineering.  ) 


There  may,  however,  be  liquid  produced  as  the  fluid  flows  up 
the  well  bore  and  into  the  separators  if  temperature  reduction 
below  the  circondentherm  occurs  within  the  tubing  string  (Path 
A-A   on  Figure  2).    If,  on  the  other  hand,  the  reservoir  is 


found  to  be  at  a  temperature  below  the  critical  temperature 
(T  )  and  pressure  above  the  bubble  point,  the  reservoir  would 
be  in  a  one  phase  liquid  state.   This  is  called  a  bubble 
point  reservoir  and  the  fluid  will  remain  a  liquid  although 
as  the  bubble  point  pressure  is  reached,  with  pressure  deple- 
tion in  the  reservoir,  gas  will  come  out  of  solution  in  the 
oil  and  will  be  produced  along  with  the  oil  in  an  ever  increas- 
ing quantities  (Path  C  -  C  ..  in  Figure  2).   In  dew  point  or  retro- 
grade gas  condensate  reservoirs  the  reservoir  fluid  exists  at 
a  temperature  between  the  critical  temperature  and  the  cricon- 
dentherm  and  at  a  pressure  above  the  dew  point  pressure  of  the 
fluid.   In  this  situation,  when  the  pressure  declines,  composi- 
tion of  the  reservoir  fluid,  which  may  be  termed  gas  in  this 
state,  will  remain  the  same  until  the  dew  point  is  reached. 
At  this  point  as  the  fluid  expands  iso therma 1 1 y ,  a  liquid  be- 
gins to  condense  out  of  the  gas  and  will  adhere  to  the  walls 
of  the  rock,  tubing,  etc.,  leaving  a  gas  leaner  in  the  heavier 
hydrocarbons  to  be  produced.   The  term  retrograde  is  given  to 
this  phenomenon  since  one  would  expect  vaporization  to  take 
place  with  a  pressure  decline  rather  than  condensation.   At 
some  intermediate  pressure,  between  the  dew  point  and  abandon- 
ment, this  retrograde  phenomenon  will  reach  a  maximum,  below 
which  pressure  the  condensate  will  revaporize  into  the  gas. 
( Pr es sure- Temp  era tur e  phase  relationships  for  gas  - c ondens a t e 
reservoirs  are  shown  on  Figure  2  by  the  path  from  B  to  B_ . ) 


B .   Retrograde  Loss  of  Condensate 

Liquid  that  is  condensed  inside  the  reservoir  rock 
remains  immobile  in  the  usual  case  and  therefore  is  lost  to 
production  unless  or  until  r evapor i za t ion  can  take  place 
either  by  further  pressure  decline,  r epr es sur i za t ion,  or  by 
sweeping  with  a  dry  gas.   The  alternative  to  r evapor i za t io n 
of  this  condensed  liquid  is  to  prevent  it  from  condensing  in 
the  first  place  and  this  possibility  leads  to  the  investiga- 
tion of  pressure  maintnenace  by  injection  of  dry  gas  into  the 
reservoir. 


C .   Experience  in  Gas  Cycling 

In  general,  the  gas  with  the  greater  condensate  content 
exhibits  the  greater  condensate  loss  upon  pressure  depletion 
so  usually  it  proves  to  be  the  more  profitable  reservoir  to 
cycle  with  gas.   In  gas  cycling  the  condensed  liquid  is  sepa- 
rated from  the  wet  (produced)  gas  in  a  gasoline  plant  or  other 
surface  facilities.   The  residue  or  dry  gas  is  then  re-injected 
into  the  reservoir  where  it  maintains  the  reservoir  pressure 
and  also  serves  to  drive  the  wet  gas  toward  the  producing  wells 
This  procedure  is  continued  until  the  dilution  by  dry  gas  at 
the  producing  well  renders  the  operation  unprofitable  at  which 
time  the  reservoir  is  blown  down,  or  allowed  to  produce  nor- 
mally through  pressure  depletion.   During  this  pressure  blow- 
down  not  only  will  the  lean  gas  in  the  reservoir  be  produced 


but  also  additional  liquid  may  be  recovered  from  both  the 
dry  gas  invaded  zone  and  the  uninvaded  portion  of  the 

res  ervo  ir . 

2 
Bowers   has  shown  that  cycling  in  the  Grapeland  Field, 

Houston  County,  Texas,  resulted  in  a  recovery  of  seventy-five 

per  cent  of  the  initial  butanes  and  heavier  hydrocarbons  in 

the  reservoir,  while  straight  pressure  depletion  would  have 

recovered  only  thirty-nine  per  cent  of  the  initial  butanes-- 

plus.   The  liquid  that  could  have  been  recovered  by  pressure 

depletion,  then,  was  only  fifty-two  per  cent  of  that  which 

was  actually  produced  by  the  cycling  process.   In  this  example 

it  was  shown  that  the  ultimate  liquid  recovery  was  reached  in 

approximately  eleven  years  but  twenty-one  years  would  have  been 

required  by  straight  depletion,  thereby  making  the  present  worth 

of  the  liquid  products  much  greater  with  gas  cycling. 

In  a  study  of  the  Krotz  Springs  Field,  St.  Landry  Parish, 

3 
Louisiana,  Robertson   predicted  that  twenty  per  cent  more  con- 
densate could  be  recovered  by  cycling,  resulting  in  a  gross 
cash  gain  after  taxes  of  $2,000,000  for  the  twenty-five  year 

period  studied. 

4 
T.  W.  Brinkley   concluded  that  benefits  in  improved 

recovery  of  stock  tank  condensate  by  pressure  maintenance 

methods  may  be  as  great  as  300  per  cent  in  rich  condensate 

reservoirs  but  grade  progressively  downward  to  negligible 

benefits  for  the  lean  condensate  reservoirs  with  GOR's 

greater  than  300,000  cubic  foot/bbl. 


Although  cycling  may  appear  to  be  the  ideal  solution 
to  the  retrograde  condensation  problem  there  are  other  fac- 
tors that  may  make  this  procedure  uneconomical.   The  deferred 
income  from  the  dry  gas  sales  may  prove  to  be  a  significant 
factor  in  the  economic  analysis  of  the  field.   Another  eco- 
nomic consideration  is  the  added  investment  required  for 
additional  injection  wells,  gas  compression  and  distribution 
system  to  return  residue  gas  to  the  wells  and  for  a  liquid 
recovery  plant.   Also,  it  must  be  remembered  that  even  with 
gas  injection  at  a  pressure  above  the  dew  point  all  of  the 
liquid  hydrocarbons  may  not  be  recovered. 

Three  different  recovery  factors  must  be  applied  with 
gas  cycling.   The  microscopic  displacement  efficiency  is 
approximately  70  to  90%.   The  volumetric  sweep  efficiency, 
or  the  per  cent  of  initial  pore  volume  invaded  by  the  sweep- 
ing gas  to  abandonment  of  the  producing  wells  can  vary  from 
50  to  9  0%.   Finally  a  permeability  stratification  factor 
must  be  applied  to  take  care  of  the  problem  of  gas  sweeping 
through  the  more  permeable  strata  and  reaching  the  producing 
wells  before  the  tighter  stringers  can  be  swept.   The  over- 
all condensate  recovery  factor  for  a  gas  cycling  project  is 
the  product  of  these  three  factors.   On  the  plus  side  of 
recovery  by  cycling  is  the  fact  that  additional  liquid  is 
produced  during  blow-down  from  the  less  permeable  or  unswept 
portions  of  the  reservoir  as  these  products  bleed  into  the 


10 

more  permeable  streaks  and  ultimately  to  the  wells.   Also  the 
efficiency  of  extracting  condensate  from  the  wet  gas  is  higher 
with  the  type  of  separator  plant  that  is  generally  used  in  a 
cycling  project. 

Any  reservoir  exhibiting  retrograde  properties  and  hav- 
ing a  richness  of  50  bbls/MMSCF  of  condensate  is  a  good  pros- 
pect for  cycling  of  the  dry  residue  gas  to  improve  recovery. 
Each  reservoir,  however,  is  an  individual  case  and,  as  such, 
requires  a  separate  study.   Variables  such  as  permeability 
stratification,  reservoir  geometry,  gas  properties,  and  market 
availability  can  team  up  to  make  one  reservoir  an  extremely 
profitable  enterprise  while  another  may  prove  entirely  infea- 
sible.   This  thesis  presents  the  study  of  one  field  based  on 
the  fluid  properties  and  geologic  structure  of  that  one  field. 


CHAPTER  III 
DEFINITION  OF  THE  FIELD  STUDIED 

A .   Development  of  the  Field 

Production  from  this  field  began  in  January,  1961  and 
from  that  time  two  to  three  wells  have  produced  continuously 
as  shown  in  Table  II.   The  pool  is  one  in  which  higher  than 
normal  pressures  are  encountered.   It  has  been  obvious  from 
the  low  gas-oil  ratios  measured  that  the  produced  gas  was 
rich  in  condensate.   By  January  1,  1964,  a  total  of  20  billion 
standard  cubic  feet  of  wet  gas  had  been  produced  from  the  two 
wells  completed  in  the  zone  studied.   The  pressure  had  dropped 
from  an  estimated  initial  pressure  of  8838  psia  to  8100  psia 
as  shown  in  Figure  7.   Based  on  these  production  data,  calcu- 
lation of  the  initial  gas  in  place  yielded  the  value  of  606.9 
billion  standard  cubic  feet.   These  values  form  the  basis  for 
future  performance  predictions  from  the  base  date  of  January, 
1964. 


B .   Structure  and  Reserves 

Interpretation  of  available  geologic  data  indicates  a 
deep-seated  salt  domal  structure.   Producing  depths  at  which 
this  study  was  conducted  were  from  approximately  11,300  to 
11,600  feet.   A  structure  contour  map  of  the  top  of  the  pro- 
ducing formation  is  presented  as  Figure  3.   The  productive 
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13 
interval  is  taken  to  be  a  shaly-sand  with  in t er gr anu  lar 
porosity. 

There  are  a  total  of  nine  wells  drilled  through  the 
formation  studied  and  from  which  data  are  available  from 
various  kinds  of  logs.   Porosity  was  determined  using  micro- 
logs  from  some  wells  and  sonic  logs  from  others.   The  porosity 
throughout  the  formation  is  fairly  constant,  varying  from  2  5 "L 
to  29%  in  the  pay  zone.   Although  shaly  streaks  were  noted, 
this  fact  was  taken  into  consideration  in  counting  sand, 
wherein  the  shale  sections  were  discounted.   For  purposes  of 
the  material  balance  calculations  an  average  uniform  porosity 
of  28 7»  was  used. 

Electric  logs  or  indue t io n- e 1 ec tr ic  logs  were  available 
from  all  nine  wells.   Three  of  these  logs  are  shown  in  Figure  4. 
These  logs  were  used  to  compute  formation  water  resistivity 
using  the  Spontaneous  Potential  Curve.   Water  saturation  was 
calculated  from  the  induction  log  or  lateral  curve  as  avail- 
able.  Water  resistivity  was  found  to  be  approximately  .03  ohm- 
meters  at  formation  temperature  and  this  value  was  used  through- 
out the  field.   Computed  water  saturation  tended  to  vary  from 
8-107o  near  the  top  of  the  gas  zone  to  about  257,  nearer  the 
bottom  of  the  pay  and  it  grades  to  10  07,  water  saturation  in 
those  wells  where  a  gas-water  contact  is  observed.   A  value  of 
207>  wa  ter- saturation  was  used  for  the  entire  formation  as  being 
the  weighted  average  value. 
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Side  wall  cores  were  available  from  one  well.   The 
porosity  values  obtained  from  the  sidewall  cores  checked  pre- 
cisely the  values  obtained  from  the  sonic  logs  and  micrologs. 
The  permeability,  however,  did  not  appear  uniform,  nor  did 
these  values  check  well  against  the  effective  permeability 
calculated  from  a  pressure  draw-down  test  or  against  permea- 
bility calculated  from  a  resistivity  gradient  observed  in  some 
of  the  wells.   The  six  values  obtained  from  side-wall  cores 
were  84,  134,  50,  512,  131  and  116  mi  1 1 i dar c i es .   The  effec- 
tive permeability  to  gas  obtained  in  the  pressure  draw-down 
test  proved  to  be  50  millidarcies  and  this  was  the  value  used 
in  the  material  balance  calculations  since  the  method  of  tak- 
ing side-wall  cores  tends  to  result  in  fractured,  non- 
representative  samples  of  the  reservoir.   A  plot  of  the  data 
taken  from  the  pressure  draw-down  test  is  shown  in  Figure  5. 
For  volumetric  sweep-out  calculations,  50  millidarcies 
was  used  for  permeability  while  the  effective  thickness  of  pay 
was  allowed  to  vary  over  the  field.   This  gave  a  value  of  flow 
capacity  for  each  point  in  the  reservoir  which  was  used  in  the 
sweep  efficiency  program. 

Connate  water  saturation  was  found  to  be  20  per  cent. 
This  was  an  average  value  based  on  water  saturation  in  the 
lower  portion  of  the  gas  zone.   Residual  gas  saturation  to 
water  is  estimated  to  be  30  per  cent.   This  is  justified  by 
the  gas  saturation  near  the  top  of  the  water  column  which,  it 


Pressure  Drawdown  Test 


Figure  5 


16 


o 
w 

<u 

S-i 
3 
co 
co 
CD 
U 

cu 

QJ 

o 

e 

o 

4J 

AJ 

o 
pq 

7900 

7875 

7850 

7825 

7800- 

7775 

7750 

7725 

\ 

! 

< 

i 
i 

! 
i 

. 

m  =   34.5 

psi/cycle 

1 

^ 

..> 

V 

1 

• 

> 

\ 

V 

Y 

v 

•01               .02        .03        .05                  .1               .2           .3           .5                    1.0 

Time    in   Days 

17 
can  be  assumed,  has  encroached  into  the  reservoir  to  some 
extent  during  a  recent  producing  time. 

Determination  of  net  effective  thickness  in  this  forma- 
tion was  very  difficult  because  of  the  shalyness  of  the  sand 
and  because  of  the  meager  geologic  information.   Only  nine 
wells  have  penetrated  the  2000  acres  of  estimated  pay  forma- 
tion.  The  locations  of  the  two  faults--one  thought  to  have  a 
closure  of  about  100  feet  in  the  northern  extremity  of  the 
pool,  the  other  with  an  estimated  closure  of  about  120  feet 
in  the  south--are  by  no  means  exactly  defined.   An  error  of 
the  position  of  the  fault  to  the  north  of  only  a  few  feet 
would  produce  a  substantial  error  in  the  estimated  reserves 
because  of  the  great  thickness  of  pay  adjacent  to  this  fault. 
Calculation  of  reserves  based  on  the  volumetric  method  is 
presented  in. Table  II.   The  east  and  west  boundaries  of  the 
field  are  located  at  a  gas-water  contact  which  has  shown 
little  or  no  encroachment  during  the  current  gas  production. 

An  isopach  map  produced  from  the  nine  control  wells  is 
presented  as  Figure  6. 

A  volumetric  gas  balance  was  made  using  the  production 
data  to  date  and  bottom  hole  pressure  tests  taken  in  1961  and 
1963.   A  bottom  hole  pressure  was  obtained  from  the  draw-down 
and  build-up  pressure  test  conducted  in  connection  with  this 
study.   This  gas  balance  was  calculated  on  a  Control  Data 
Corporation  1604  high  speed  computer  using  a  program  developed 
by  the  Union  Oil  Company  of  California. 
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J.  Donald  Clark,  Chief  Reservoir  Engineer  for  the  Gulf 
Division  of  the  Union  Oil  Company  of  California,  has  developed 
and  applied  or  supervised  the  application  of  a  very  simple 

transient  energy  balance  equation  to  numerous  natural  water 

12 
drive  gas  reservoirs  found  in  the  Gulf  Coast  area.    A  primary 

requirement  of  this  equation  is  the  derivation  of  the  gas  for- 
mation volume  factor,  B  ,  to  convert  standard  cubic  feet  of 
gas  to  reservoir  barrels.   This  conversion  for  any  particular 
reservoir  is  presented  mathematically  as  follows: 


B   ■  15-°25  Tf  x  |  -  .0051505  Tf  Z/P 
8    520(5.61) 


Where:   15.025  is  base  pressure  for  Louisiana  (psi) 
520  is  standard  temperature  °R 
5.61  converts  cubic  feet  to  barrels 
T  is  reservoir  temperature,  degrees  Rankine 
Z  is  gas  deviation  or  compressibility  factor 
P  is  the  reservoir  pressure,  psia 
i  as  a  subscript  will  denote  initial  or 
original  conditions 


The  cumulative  reservoir  barrels  of  gas  produced  or  withdrawn 
is  equal  to  the  cumulative  reservoir  barrels  of  gas  expansion 
in  the  reservoir,  or 


G    B 

pw   g 


Wher  e 


G(B  -B  .  ) 

g   gi 


Gpw  is  cumulative  standard  cubic  feet  of 

wet  gas  produced 
G  is  standard  cubic  feet  of  original  gas 

in  place 


In  the  case  of  water  drive  reservoirs  and  water  production 
the  statement  of  the  equation  then  becomes:    (Cumulative  Net 
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Reservoir  Withdrawals)  -  (Cumulative  Reservoir  Expansion)  = 
(Cumulative  Water  Influx),  or 

3.  (G    B   +W)-G(B   -  B  .  )  =  W 

pw   g     p         g     gi      e 

4.  (G    B   +  W  )  -  G(B   -  B  .)  =  AEAPQ 

pw   g      p         g     gi  r 

The  water  influx  equation,  W   =  AEAPQ  ,  is  a  simplified 
method  of  writing  the  equation  presented  by  Hurst   in  his 
paper  entitled  "Water  Influx  into  a  Reservoir  and  Its  Appli- 
cation to  the  Equation  of  Volumetric  Balance."   Details  of 
solving  the  water  influx  equation  were  presented  by  van 
Everdingen  and  Hurst.    The  "A"  is  the  water  influx  constant 
which  encompasses  conversion  factors,  porosity,  radius  of 
reservoir,  effective  reservoir  thickness,  effective  compressi' 
bility  of  reservoir  fluid,  viscosity  of  reservoir  fluid  and 
fraction  of  reservoir  perimeter  exposed  to  water  drive.   The 
above  parameters  are  constant  for  any  particular  reservoir 
and  when  used  with  the  volumetric  balance  equation  can  be 
solved  as  an  all  encompassing  constant  for  that  particular 
res  ervo  ir . 

The  true  simplification  is  in  the  manner  of  arriving  at 

the  value  of  EAPQ   for  various  balancing  periods. 

2 

The  dimens ionl es s  time  equation,  t^  =  KT/0mC  r,    is 

n        '   D  w  b 

2 
derived  from  the  hydraulic  diffusivity  factor,  a   =  K./pa'C 

13 
presented  by  Hurst   in  1933.   The  K,  or  permeability  factor, 

actually  is  the  one  most  difficult  to  measure  since  it  is 


21 
effective  permeability  to  water,  which  we  seldom  if  ever 
obtain.   This  value  could  be  found  by  running  draw-down  tests 
in  salt  water  wells,  but  actually  it  can  be  estimated  to  a 
sufficient  degree  of  accuracy.   Viscosity  of  water  can  be 
obtained  from  various  sources  to  a  reasonable  order  of  magni- 
tude.  Therefore,  it  becomes  important  only  to  establish  an 
order  of  magnitude  value  of  t  ,  in  order  to  extract  reasonable 
values  of  the  corresponding  Q  ,  the  dimens ion  1 es s  water  influx, 
from  the  tables. 

Experience  has  shown  that  a  proper  time  period  should 
not  exceed  three  months  or  approximately  a  90-day  producing 
period  when  attempting  to  strike  a  balance,  especially  if  one 
plans  to  calculate  future  predictions.   The  constant  "A"  be- 
comes an  over-all  correcting  constant  of  transient  flow,  and 
allows  us  to  solve  very  simply  a  satisfactory  value  of  original 
hydrocarbons  in-place  and  the  total  water  influx  at  any  time  in 
the  producing  life  of  the  reservoir. 


Cumulative  Production  and  Pressure  History 

Production  from  the  reservoir  is  tabulated  on  a  monthly 
rate  and  cumulated  on  a  three-month  basis.    These  data  are 
broken  down  between  oil,  gas,  condensate  and  water  (both 
fresh  and  salt  water).   When  rates  of  production  have  consider 
able  variation,  a  monthly  production  rate  should  be  plotted  as 
shown  on  Figure  7.   This  plot  is  important  in  order  to  show 
the  effect  of  rate  on  pressure  when  natural  water  drives  exist 
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Bottom  hole  pressure  measurements  are  plotted  versus 
time  on  the  same  graph  as  production  rate.   The  length  of 
shut-in  time  prior  to  each  pressure  measurement  should  be 
approximately  equal  and  the  individual  pressures  should  have 
reasonably  close  values.   If  not,  build-up  or  drawdown  tests 
should  be  run  to  determine  the  cause  of  pressure  variation. 
A  pressure  curve  is  then  constructed  through  these  evaluated 
pressure  points.   Experience  and  ingenuity  are  of  tremendous 
value  in  constructing  a  reasonable  curve.   The  production 
rate  graph  must  be  studied  in  order  to  prepare  the  most  logi- 
cal pr es sur e- t ime  curve.   Pressure  for  uniform  periods  of  pro- 
duction are  taken  directly  from  this  graph. 

Water  Influx  Data  and  Basic  Calculations 

The  next  step  requires  the  gathering  of  basic  informa- 
tion to  approximate  d imens ion  1 es s  time  data  and  corresponding 
dimens ion  1 es s  rates. 


KT(86,400) 


D     u0(O. 5  x  10"4)  r, 2(930) 


Here  geological  data  can  be  used  to  help  determine  the  radius 
of  the  hydrocarbon  reservoir.   Only  a  close  approximation  is 
necessary.   Compressibility  of  salt  water  can  be  found  in  the 
literature.   Viscosity  of  salt  water  can  be  approximated  by 
using  chemical  handbooks,  available  texts,  and  various  publi- 
cations; porosity  is  approximated  from  core  analyses  or  logs 
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and  permeability  to  water  can  be  estimated  from  air  permea- 
bility of  cores  analyzed. 

Dimens ionl es s  time  values  are  calculated  for  cumulative 
three-month  periods.   The  corresponding  Q   values  are  taken 

from  the  table.   Integration  of  the  values  for  AP  times  Q 

r 

over  the  cumulative  time  intervals  yields  the  value  which, 
when  multiplied  by  the  constant  A;  provides  the  cumulative 
wa  t  er  inf  lux . 

Gas  Conversion  to  Reservoir  Barrels 

The  vapor  equivalent  of  the  distillate  and  fresh  water 
produced  must  be  added  to  the  amount  of  dry  gas  produced.   The 
constant  for  converting  liquid  water  to  the  vapor  state  is 
7368.48  cubic  feet  per  barrel.   The  distillate  conversion  fac- 
tor was  calculated  from  the  hydrocarbon  analysis  of  the  sepa- 
rator liquid  shown  on  page  89,  Appendix  II.   This  factor  was 
found  to  be  768  standard  cubic  feet  of  vapor  per  barrel  of 
distillate. 

Method  of  Least  Mean  Squares  for  Striking  an 
Energy  Ba lance 

The  method  of  least  mean  squares  is  used  for  striking 

an  energy  balance  for  the  natural  water  drive  gas  reservoir. 


(G    +W   )  -  G(B   -B.)=  A(EAPQ  ) 
pw     ps;        g     giy  xr 

The  two  unknowns  are  G,  the  original  wet  gas  in-place, 
and  A,  the  water  influx  constant.   The  method  herein  presented 


25 

requires  the  solution  for  single  values  of  A  and  G  to  meet  all 
periods  of  past  producing  time.   We,  therefore,  can  arrive  at 
two  equations  with  two  unknowns  and  place  proportional  weight 
to  the  length  of  production  time: 

[E(G    +  W   )(EAPQ  )]  -  G[(B   -  B   ) (EAPQ  ) ]  =  A[(EAPQ  )2] 

[E(G    +  W   )(B   -  B  .)]  -  G[(B   -B.)2]=  A  [(EAPQ  )(B  -B  .)] 
pw     ps    g     gi/J        g     gi'J  xr'v  g   gi/J 

These  two  equations  can  then  be  solved  for  A,  the  water 
influx  constant,  and  G,  the  original  wet  gas  in-place. 

In  order  to  check,  the  validity  of  the  values  of  A  and 
G  obtained  by  this  method,  multiply  the  LMSA  value  of  A  times 
EAPQ   and  solve  for  the  apparent  value  of  G.   This  is  called 
"periodic  check  of  balance  for  G  assuming  LMSA  values  of  A 
as  constant."   The  van  Everdingen  and  Hurst  solution  to  the 
diffusivity  equation  is  here  applied  to  single-phase  fluid 
transient  flow  into  a  well  bore  for  determination  of  the  rate 
and  volume  of  encroachment  of  water.   The  gas  field  is  assumed 
to  be  embedded  into  an  aquifer  of  large  extent.   In  this 
instance,  the  reservoir  is  taken  to  be  a  large  well,  the 
radius  of  which  is  the  equivalent  field  radius,  and  the  aqui- 
fer is  considered  to  be  the  reservoir  from  which  fluid  is 

g 
withdrawn  and  water  is  expelled. 

The  computer  program,  indicating  method  of  computation, 

reservoir  conditions,  initial  gas  in  place,  and  cumulative 


26 

water  influx  from  past  production  data,  is  included  as 
App  end  ix  I . 

Since  the  available  geologic  information  is  considered 
isufficient  to  describe  the  reservoir  accurately  enough  for 
calculation  of  reserves,  the  results  obtained  from  the  volu- 
metric gas  balance  were  used  for  estimating  original  gas  in 
place  and  water  encroachment.   These  results  were  also  the 
basis  for  the  future  performance  prediction  and  economic 
analysis  of  alternate  methods  of  exploitation.   The  original 
gas  in  place  was  found  to  be  606.9  billion  cubic  feet  measured 
at  standard  conditions  (15.025  psi  at  60  F).   Water  influx  was 
found  to  be  negligible. 


C .   Reservoir  Fluid  Characteristics 

Samples  of  both  the  separator  liquid  and  gas  were  ob- 
tained in  March,  1964,  and  recombined  in  accordance  with  the 

measured  gas-liquid  ratio  (stock-tank)  of  17,150  SCF  separator 

«. 
gas  per  bbl.  stock  tank  liquid.    It  was  difficult  to  obtain 

constant  gas-oil  ratios  for  the  purpose  of  getting  a  repre- 
sentative sample  and  flowing  at  a  constant  rate  for  about  three 
days  was  required  before  a  steady  GOR  was  achieved.   Analysis 
of  the  recombined  sample  indicates  a  dew  point  pressure  above 
the  present  bottom  hole  pressure.   Retrograde  condensation 
was  taking  place  and  two  phase  flow  occurring  near  the  well 
bore  caused  the  fluctuating  producing  gas-liquid  ratio.   Reser- 
voir temperature  at  time  of  sampling  was  240  degrees  farenheit. 
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Primary  separation  temperature  was  125  degrees  farenheit  and 
primary  separator  pressure  was  1020  psig.  Dew  point  pressure 
was  found  to  be  8260  psia  and  gas  richness  was  66.07  bbls.  of 
stock  tank  liquid  per  million  standard  cubic  feet  of  well 
stream  effluent.  Detailed  results  of  the  P-V-T  investigation 
are  included  as  Appendix  II. 

Samples  of  separator  gas  and  separator  liquid  were  ana- 
lyzed by  the  conventional  low  temperature  distillation  method. 
The  results  are  reported  showing  both  the  composition  of  each 
sample  and  the  computed  analysis  of  the  well  stream  based  on 
the  gas-liquid  ratio  of  the  primary  separator.  The  separator 
liquid  production  was  calculated  from  the  measured  stock  tank 
production  by  applying  the  determined  shrinkage  factor. 

Following  the  compositional  analyses,  portions  of  the 
primary  separator  liquid  and  gas  were  physically  recombined 
in  their  produced  ratio  in  a  variable  volume,  g  1  as s -windowed 
equilibrium  cell.   Determinations  on  this  mixture  were  divided 
into  the  following  two  main  categories: 

I.   Dew  point  pressure  determination  and  pressure- 
volume  relations  on  a  constant  weight  of  reservoir  fluid 
at  the  reservoir  temperature:   The  procedure  consisted  of 
establishing  equilibrium  between  gas  and  liquid  phases  at 
a  low  pressure  and  measuring  the  volumes  of  liquid  and 
gas  in  equilibrium  at  that  pressure.   The  pressure  was 
then  raised  by  the  injection  of  mercury  into  the  cell  and 
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phase  equilibrium   established  again  at  the  higher  pressure. 
This  procedure  was  repeated  until  all  of  the  liquid  phase 
had  vaporized,  at  which  point  the  saturation  pressure  was 
observed.   The  cell  pressure  was  then  raised  above  the  dew 
point  pressure  in  order  to  determine  the  sup er compr es s i- 
bility  characteristics  of  the  single  phase  vapor.   As  a 
check,  on  all  readings  and  particularly  to  verify  the  dew 
point,  the  cell  pressure  was  incrementally  reduced,  equi- 
librium established  and  volumetric  readings  made.   Reported 
in  Appendix  II  are  the  relative  volume  relations  and  spe- 
cific volumes  of  the  reservoir  fluid  over  a  wide  range  of 
pressures  as  well  as  deviation  factors  of  the  single  phase 
vapor  above  the  dew  point.   (Also  reported  in  Appendix  II 
are  the  dew  point  pressures  resulting  from  recombinations 
at  gas-liquid  ratios  above  and  below  the  ratio  measured 
at  the  time  of  sampling.) 

II.   Compositions  of  the  produced  well  stream  and  the 
amount  of  retrograde  condensation  resulting  from  a  stepwise 
differential  depletion:   This  procedure  consisted  of  a 
series  of  constant  composition  expansions  and  constant  vol- 
ume displacements  with  each  displacement  being  terminated 
at  the  original  cell  volume.   The  gas  removed  during  the 
constant  pressure  displacement  was  charged  to  the  low 
temperature  fractional  distillation  equipment  for  analysis. 
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D  .   Conf ormanc  e 

The  description  of  the  simultaneous  flow  of  fluids 
through  porous  media  in  terms  of  relative  permeability  and 

capillary  pressure  has  been  adequately  discussed  in  the  lit- 

9 
erature.   For  example,  Muskat   in  Chapter  VII  briefly  dis- 
cussed the  dynamical  effects  associated  with  capillary 
phenomena.   He  also  outlined  the  theory  of  po t en t iome tr ic 
models  and  illustrated  the  application  of  the  po t en tiome tr i c 
model  to  cycling  systems.   Douglas,  Peaceman,  and  Rachford, 
presented  a  method  for  calculating  multi-dimensional  immiscible 
displacement  based  on  a  numerical  solution  of  a  finite  differ- 
ence analogue  describing  the  displacement  process.   This  tech- 
nique insures  that  at  any  point  in  the  sand,  the  flow  of  oil 
per  unit  area  across  the  direction  of  flow  can  be  represented 
by  a  vector  u   which  is  given  by 
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u   =  -K  - 

°         u 


ro 


(VpQ  +  PQ  g  Vh). 


The  operator,  v>  equals   i  -v—  +  j  -*r—  ;  i,  J  are  unit 
vectors  in  the  x  and  y  directions,  respectively;  k  is  the 
local  permeability,  k    is  the  relative  permeability  to  oil, 
u   the  oil  viscosity,  p   the  pressure  in  the  oil  phase,  p   the 
oil  density,  and  g  the  acceleration  of  gravity.   Similarly, 
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Fluid  densities  are  assumed  constant  and  oil  and  water 
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potentials  are  defined 


|   =  p   +  p   gh 

(5   =  p   +  p   gh 
xw     w     w 

The  author  proceeds  to  solve  the  problem  by  a  system  of  dif- 
ference equations. 

In  solving  this  reservoir  a  technique  demonstrated  by 
Dr.  S.  J.  Pirson  at  The  University  of  Texas  is  followed.   In 
this  method  the  reservoir  is  divided  into  a  grid  and  each 
intersection  of  the  grid  is  treated  as  a  point  (I,  J)  within 
the  reservoir  where  Darcy's  law  may  be  applied.   The  flow 
across  the  grid  between  these  points,  q,  is  now  analyzed. 
The  rate  of  flow  over  a  distince  Ax  across  an  area  open  to 
flow,  Ay,  can  be  determined  by  Darcy's  law: 
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In  order  to  satisfy  the  continuity  equation  the  vector 
sum  of  the  flow  in  all  four  directions  must  equal  zero  unless 
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there  is  fluid  withdrawal  or  fluid  injection  at  the  point 


Q  =  qx  +  q2  +  q3  +  q4 

A    =    1;127     [K.         (P.     .      T-P.      .  )  +  K,         (P.      ,^,-P,      .  )  +  K,         (P...      .-P.      .) 
0u  h  i-l. J       i,j         h  i,J+l       ijj         h  i+l, J       iij 
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It  is  seen  that  without  injection  or  production  from  the 

point  studied  0=0  and  the  Kh    is  the  average  flow  potential 

11  a  v 


between  adjacent  points.   This  equation  is  now  solved  for  P 


i>  J 


and  we  find: 
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The  boundary  conditions  are  those  restrictions  on  the 
limits  of  the  connected  region  and  the  position  and  magnitude 
of  the  imposed  flows  or  pressures.   At  the  boundary  of  the 
permeable  sand,  the  velocity  of  flow  perpendicular  to  the 
boundary  is  zero,  so  it  is  sufficient  to  establish  equal  poten- 
tials across  the  boundary.   Since  (V$)  =  0  across  the  boundary, 
no  flow  occurs.   This  same  condition  occurs  at  a  fault  bound- 
ary.  As  the  water  saturation  increases  near  the  water  table, 
the  relative  permeability  to  water  increases  greatly  while 
that  to  oil  diminishes.   By  assigning  these  relatively  high 
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values  of  flow  potential  in  the  equation  aquifer  pressure  is 
maintained  in  the  advancing  flood  front  and  movement  of  the 
front  is  regulated  across  this  boundary. 

In  solving  this  potential  equation  for  an  oil  reservoir, 
the  term   -  2|3  u  / 1  .  12  7   is  treated  as  a  constant.   However, 
since  the  flow  equation  depends  on  pressure,  the  constant  for 
a  gas  field  is  augmented  by  a  variable  pressure  term: 


q  =  - 


1.127  ^ 

|_i  dx 


For  Oil 


(q   )(P   )(T)(z) 

s  c    s  c 

(5.615)(Tsc)(A)(P) 


For  Gas 


Volumetric  sweep  efficiency  was  investigated  by  this 
technique  with  the  aid  of  a  mathematical  model  which  was  run 
on  the  CDC  1604  high  speed  computer.   This  computer  program 
is  included  as  Appendix  III.   The  program  assumes  steady-state 
flow  at  a  point,  and  sequentially  at  each  point  in  the  reser- 
voir in  order  that  Darcy's  flow  equation  may  be  applied.   Using 
an  alternating  direction  iteration  procedure,  the  computer  then 
reduces  the  total  flow  at  each  point  in  four  directions  to  zero 
and  by  many  iterations  establishes  each  potential  from  the 
adjacent  four  potentials  in  order  to  achieve  a  balance.   At 
each  point  selected  on  the  model  the  quantity  of  flow  into  or 
out  of  the  reservoir  is  assigned  in  the  event  a  well  is  located 
at  that  point.   The  flow  potential,  or  the  product  of  effective 
permeability  times  the  thickness  of  effective  pay,  is  also 
assigned  to  each  point.   By  applying  Darcy's  law  in  each  of 
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four  directions,  by  summation  of  these  directional  flow  rates 
to  zero,  and  by  solving  this  equation  for  the  potential  at 
each  point,  a  potential  distribution  map  of  the  reservoir   is 
produced.   Isopotential  contour  lines  are  now  drawn  on  a  scale 
map  of  the  reservoir.   The  primary  criterion  for  equivalence 
between  the  model  and  the  actual  reservoir  is  the  creation  of 
geometrically  similar  potential  fields  which  are  determined 
only  by  the  thickness,  permeability,  and  boundary  conditions. 
Flow  lines  were  superimposed  on  the  potential  map  and  the 
method  of  squares   utilized  by  Seba    is  used  to  define  the 
flow  tubes  through  which  a  constant  quantity  of  reservoir 
fluid  would  be  expected  to  flow. 

Flow  lines  of  a  fluid  in  a  porous  medium,  as  defined 
in  Darcy's  law,  are  everywhere  orthogonal  to  the  isopotential 
lines.   The  method  of  squares  employs  the  fact  that  if  a  family 
of  curves  is  constructed  such  that  they  are  orthogonal  to  one 
member  of  another  family  of  curves,  then  the  second  family  of 
curves  may  be  drawn  orthogonal  to  the  first.   This  is  done  by 
constructing  them  such  that  adjacent  distances  between  the 
streamlines  are  equal  to  adjacent  distances  between  the  iso- 
potentials.   A  stream  tube  is  a  two-dimensional  flow  system 
defined  as  the  area  between  two  consecutive  stream  lines.   Flow 
rate  can  be  defined  by  the  following: 


Q  =  (h)(-  -)(Ax)(^J) 

u       As 
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In  any  one  square  Ax  =  As  and  Q  is  constant  so 


(As) 


At  =  (constant)  -*-  -^ 


The  increment  of  time  can  now  be  calculated  for  each  of  the 
squares  within  each  flow  tube  so  that  the  position  of  the 
front  in  each  flow  tube  can  be  drawn  at  any  selected  time 
period.   This  procedure  is  repeated  for  various  injection- 
production  well  combinations  until  the  maximum  volume  of  the 
reservoir  is  swept.   Rearrangement  of  the  well  positions  in 
the  reservoir  matrix  is  achieved  simply  by  changing  the  quan- 
tity of  flow  into  or  out  of  that  point  in  the  reservoir.   Sweep 
out  patterns  were  predicted  with  one,  two  and  three  injection 
wells  and  by  using  four  producing  wells.   The  locations  of 
existing  wells  were  used  as  much  as  possible  but  alternate 
new  locations  were  also  tried.   The  optimum  combination  was 
found  to  be  two  injection  wells  and  four  producing  wells. 

Wells  marked  "Well  4"  and  "Well  10"  on  Figure  8  are  to 
be  recompleted  as  injection  wells.   Two  new  wells  are  required 
at  the  locations  marked  "Well  5"  and  "Well  11."   An  existing 
well  at  location  "Well  9"  is  to  be  recompleted  in  the  pay  zone 
studied  since  it  is  now  producing  from  a  deeper  zone.   "Well  6" 
is  producing  from  this  formation  and  will  continue  to  be  used 
for  this  project. 

A  machine  print-out  of  the  potentials  has  been  reduced 
to  a  map  shown  in  Figure  8  picturing  the  isopotential  contour 
lines  and  on  which  the  frontal  position  has  been  drawn  at  one 
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time  just  after  breakthrough  of  the  injected  residue  gas.   The 
volume  of  the  reservoir  swept  at  this  breakthrough  time  was 
found  to  be  65%. 

The  ability  of  the  dry  absorber  gas  to  vaporize  the 
condensate  in  the  sand  was  determined  by  pr es sur e- vo lume- 
temperature  tests  run  on  a  recombined  sample  of  the  gas. 
Results  of  these  tests  are  given  in  Appendix  II.   The  overall 
displacement  efficiency  of  the  residue  gas  was  found  to  be 
9  2%. 

The  effect  of  bypassing  on  a  local  scale,  due  to  permea- 
bility variations,  was  analyzed  as  demonstrated  by  Standing, 
by  arriving  at  a  permeability  variation  for  the  sand.   Although 
insufficient  data  are  available  to  establish  firmly  the  perme- 
ability variations  throughout  the  reservoir,  a  variation  .30 
was  assigned.   This  value  gives  a  sweep  efficiency  at  break- 
through of  87%. 

Using  the  product  of  these  three  conformance  factors, 
the  overall  volume  of  the  reservoir  swept  out  at  breakthrough 
was  51%,  with  the  dry  gas  cut  after  breakthrough  increasing  in 
the  manner  calculated  by  Standing   in  the  paper  mentioned  before. 
Cycling  of  gas  was  discontinued  when  proportion  of  reservoir  gas 
produced  was  reduced  to  677«  by  dilution  with  cycled  residue  gas, 
due  to  the  cycling  economic  limit  having  been  reached. 


CHAPTER  IV 
PERFORMANCE  PREDICTION 

A.    Pressure  Depletion  Method 

Past  production  has  averaged  25  million  standard  cubic 
feet  per  day  from  the  reservoir  studied.   Until  this  time  no 
firm  estimate  had  been  made  on  the  reserves  nor  had  any  analy- 
sis been  made  to  predict  the  future  performance  and  probable 
life  of  the  reservoir.   From  the  past  production  figures  the 
initial  gas  in  place  was  estimated  and  the  pressure  decline 
and  water  influx  was  analyzed.   From  these  data  it  is  possible 
to  predict  the  future  performance  of  the  field.   This  was  done 
through  the  use  of  the  computer  program  included  in  Appendix  IV 

Because  of  insufficient  pressure  decline  and  the  lack  of 
back  pressure  data  it  is  not  possible  to  predict  de  1  i ver ab i  1  i ty 
deterioration  with  declining  pressure.   It  is  anticipated  that 
as  liquid  saturation  around  the  well  bore  increases  with  static 
pressure  reduction  and  with  more  severe  pressure  drawdown, 
deterioration  of  de 1 iver abi 1 i ty  will  result.   For  the  purposes 
of  this  study,  in  absence  of  definitive  data,  the  rate  of  pro- 
duction of  dry  gas  was  assumed  to  remain  constant  and  at  the 
current  rate.   The  water  influx  constant  obtained  from  the 
volumetric  gas  balance  program  was  used  to  obtain  future  per- 
formance prediction  as  well  as  the  calculated  quantity  of 
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initial  gas  in  place  and  the  adjusted  pressure  decline 
figures . 

The  amount  of  condensate  produced,  by  year,  was  calcu- 
lated from  the  P-V-T  data  obtained  in  the  material  balance 
calculations  performed  on  the  recombined  gas  sample  placed 
in  the  equilibrium  test  cell.   In  this  cell  the  pressure  was 
depleted  by  removing  successively  small  quantities  of  gas  from 
the  cell.   Thermodynamic  equilibrium  was  established,  material 
balance  calculations  were  made,  and  retrograde  liquid  volumes 
and  the  gravity  of  the  effluent  gas  were  measured  at  each 
step.   The  gas  removed  in  each  depletion  step  was  charged  to 
the  low  temperature  distillation  equipment  for  analysis.   The 
produced  volumes  of  reservoir  gas,  the  determined  compositions, 
computed  gallons  per  thousand  standard  cubic  feet  (GPM)  con- 
tent and  respective  deviation  factors  are  presented  in  Appen- 
dix II.   The  volume  of  retrograde  liquid  resulting  from  this 
depletion  study  is  also  shown  in  Appendix  II,  both  in  terms 
of  barrels  of  reservoir  liquid  and  per  cent  of  hydrocarbon 
pore  space.   Shown  also  are  the  compositions  of  the  gas  and 
liquid  remaining  in  the  reservoir  after  depletion  to  abandon- 
ment pressure.   The  amount  of  condensate  calculated  to  be  pro- 
duced with  the  well  stream  at  each  pressure  decrement  in  the 
future  prediction  was  tabulated.   The  value  of  gas  sold  under 
the  terms  of  the  present  contract  is  $.20  per  MC F  and  conden- 
sate is  sold  for  $2.95  per  barrel.   These  prices  are  assumed 
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to  be  firm  for  future  sales.   Operating  costs  are  averaging 
$.00423  per  MC F  produced  and  product  transportation  charges 
for  the  condensate  sold  are  averaging  $.1375  per  barrel. 
Application  of  these  values  to  depletion  of  the  reservoir 
over  a  fifty-year  period  with  net  income  discounted  at  6  % 
results  in  a  present  worth  of  $38,749,249.   It  is  seen  that 
in  this  case  of  pressure  depletion  to  1500  psi  only  73.51  per 
cent  of  the  initial  gas  in  place  was  produced.   Since  the  gas 
is  to  be  delivered  to  a  pipeline  at  a  pressure  in  excess  of 
1000  psi,  further  recovery  would  require  recompression  to  the 
sales  gas  line.   This  secondary  recovery  process  was  not  con- 
sidered in  this  study  because  the  value  of  gas  recovered  50-odd 
years  hence  is  not  relevant  for  a  study  of  possible  investments 
now . 

It  is  significant  to  note  that  the  condensate  recovered 
by  pressure  depletion  amounts  to  15.4  million  barrels  which 
is  scarcely  38%  of  the  more  than  40  million  barrels  of  recov- 
erable oil  initially  in  place. 

B ,   Pressure  Maintenance  through  Gas  Cycling 


A  program  of  gas  cycling  was  considered  wherein  gas 
sales  from  the  reservoir  would  be  terminated  and  all  gas  pro- 
duced from  the  reservoir,  together  with  enough  make-up  gas  to 
allow  for  compression  and  other  gas  losses,  would  be  reinjected 
after  extraction  of  all  condensable  heavier  hydrocarbons.   As 
determined  by  the  sweep  efficiency  model  studies,  two  wells 
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should  be  recompleted  as  injection  wells  and  four  wells  would 
be  used  as  producing  wells  as  shown  in  Figure  8.   Since  the 
installation  is  located  off  shore  the  wet  gas  would  be  pro- 
cessed in  low  temperature  units  on  platforms  near  the  produc- 
ing wells.   The  residue  gas  would  be  piped  from  the  separators 
to  the  injection  platform  which  would  be  located  at  one  of  the 
injection  wells.   Double  extra  heavy  pipe  would  be  necessary 
to  carry  the  compressed  injection  gas  to  the  other  injection 
wells  since  injection  pressures  in  excess  of  8000  psi  will  be 
encoun  t  er  ed . 

Condensate  produced  in  the  separator  units  would  be 
metered  and  then  reinjected  into  the  gas  sales  line  of  the 
gas  company  where  it  would  be  transported  with  other  wet  gas 
to  the  on-shore  extraction  plant  for  processing. 

Under  pressure  maintenance  it  is  assumed  that  the  reser- 
voir will  continue  to  produce  a  gas  of  the  richness  measured 
in  the  most  recent  flow  tests  until  breakthrough  of  the  in- 
jected dry  gas  into  the  first  producing  well.   From  this  point 
in  time  the  produced  gas  would  decline  in  richness  in  direct 
proportion  to  the  residue  gas  cut  produced.   When  the  amount 
of  injected  dry  gas  obtained  through  the  production  wells 
reached  1/3  of  the  total  gas  processed  it  was  found  uneconomical 
to  continue  gas  injection.   Injection  would  be  discontinued  and 
gas  sales  would  begin  to  deplete  the  reservoir  to  about 
1300  psi  or  until  compression  was  necessary  to  charge  the  sales 
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gas  line.  Production  past  this  point  was  not  considered  part 
of  this  study. 

It  was  assumed  that  the  wet  gas  produced  would  be  in 
equilibrium  at  the  pressure  and  temperature  produced  in  the 
same  proportion  as  that  indicated  by  the  equilibrium  cell 
tests.   In  the  equilibrium  cell  the  recombined  gas  sample  was 
allowed  to  expand  at  constant  composition  with  each  pressure 
decrement  and  then  gas  was  extracted  at  constant  pressure, 
thus  simulating  reservoir  withdrawals.   The  material  balance 
calculations  resulting  from  these  tests  are  believed  to  depict 
more  realistically  the  actual  reservoir  performance  than  calcu- 
lation of  the  material  balance  by  means  of  equilibrium  constants 

Based  on  a  permeability  variation  of  .30  as  defined  by 
Standing   the  dry  gas  cut  at  the  producing  wells  was  estimated 
from  the  first  breakthrough  at  a  producing  well  until  the  eco- 
nomic limit  of  injection  was  reached. 

Extraneous  gas  is  available  in  excess  of  present  avail- 
able market  in  a  reservoir  below  the  zone  studied  and  conse- 
quantly  at  a  higher  pressure.   This  gas  would  be  used  as  make 
up  gas  in  the  pressure  maintenance  case  to  compensate  for 
losses  due  to  fuel  consumption  and  condensate  production. 

After  the  economic  limit  is  reached  it  is  assumed  that 
gas  sales  will  begin  again  at  the  present  rate  and  that  pre- 
vailing prices  would  still  be  applicable. 

Since  the  field  is  just  beginning  to  enter  the  retro- 
grade condensation  phase  considerably  more  condensate  will  be 
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produced  by  pressure  maintenance  if  this  procedure  can  be 
started  at  once.   The  greater  the  decline  in  pressure  before 
injection  commences  the  greater  will  be  the  losses  in  the 
reservoir  and  the  higher  will  be  the  value  of  extraneous  gas 
inj  ec  t ed . 

Because  of  the  loss  in  gas  sales  during  the  period  of 
pressure  maintenance  this  procedure  does  not  appear  economi- 
cally attractive  unless  market  availability  for  present  gas 
sales  can  be  filled  from  another  source  of  gas.   The  ultimate 
recovery  of  condensate  through  pressure  maintenance  is  excel- 
lent from  a  conservation  point  of  view  since  over  907o  of  the 
available  condensable  hydrocarbon  is  expected  to  be  recovered 
through  the  combination  of  dry  gas  injection  to  maintain  pres- 
sure,  to  be  followed  by  reservoir  blowdown  to  depletion. 


C .   Cycling  with  Pressure  Depletion 

Return  to  the  reservoir  of  available  gas  from  the  high 
pressure  separator  constitutes  normal  cycling.   Since  a  net 
voidage  situation  exists  as  a  result  of  shrinkage  and  fuel 
losses  from  the  produced  gas,  normal  cycling  results  in  only 
partial  pressure  maintenance.   However,  the  pressure  reduc- 
tion would  be  little  more  than  407o  during  the  full  economic 
cycling  life. 

As  in  the  case  studied  for  full  pressure  maintenance, 
two  wells  would  be  used  for  injection  of  residue  gas  and  four 
wells  would  be  used  for  production  of  the  wet  gas  effluent  as 
indicated  in  the  sweep  model  study.    (See  Figure  8.) 
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In  this  analysis  gas  sales  are  continued  at  the  present 
rate  which  is  dictated  by  the  available  gas  market.   Provision 
is  made  for  compressors  adequate  to  compress  40  million  stand- 
ard cubic  feet  of  residue  gas  per  day  for  injection.   All  wet 
gas  produced  is  processed  through  the  low  temperature  units  on 
offshore  platforms.   The  sales  gas  and  all  the  condensate  pro- 
duced is  charged  to  the  sales  gas  line  provided  by  the  gas 
company . 

The  amount  of  condensate  produced  will  be  governed  by 
the  equilibrium  saturation  at  the  prevailing  producing  bottom 
hole  pressure  until  breakthrough  of  detectable  dry  gas.   From 
this  point  onward  the  increasing  dry  gas  cut  will  be  governed 
by  the  permeability  stratification  as  well  as  by  the  areal 
sweep  pattern.   As  the  cycling  reaches  its  economic  limit,  the 
compressors  will  be  shut  down  and  sales  gas  production  will 
continue  to  deplete  the  field.   Production  below  1300  psi  was 
not  considered  in  this  study,  since  compression  would  be  re- 
quired to  charge  the  sales  gas  line  which  is  maintained  in 
excess  of  1000  psi. 

It  is  predicted  that  normal  cycling  to  .67  reservoir- 
gas  cut  followed  by  pressure  blowdown  will  result  in  the  recov- 
ery of  26,876,000  stock  tank  barrels  of  condensate  or  67%  of 
the  available  condensable  hydrocarbons.   Although  the  total 
amount  of  condensate  produced  through  normal  cycling  is  con- 
siderably less  than  for  the  pressure  maintenance  case,  the 
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present  worth  of  the  income  from  early  gas  sales  tended  to 
make  normal  cycling  more  attractive  economically.   The  ad- 
vantage of  normal  gas  cycling,  then,  is  continued  income  from 
the  available  market  for  gas  while  producing  more  condensate 
than  would  have  been  produced  through  normal  pressure  deple- 
tion and  by  producing  this  condensate  earlier  in  the  life  of 
the  field  and  in  excess  of  available  market  for  residue  gas. 


CHAPTER  V 
CONCLUSIONS 


A .   Comparison  of  Ultimate  Recoveries 

In  a  condensate  reservoir,  income  is  derived  from  the 
sales,  not  only  of  the  residual  gas  produced,  but  also  the 
liquids  produced  from  the  reservoir  gas  by  mechanical  sepa- 
ration or  other  type  of  plant  processing.   The  ultimate 
quantity  of  these  fluids  produced  is  of  importance  as  well 
as  the  rate  of  this  production. 

The  original  gas  in  place  is  the  same  in  any  case. 
The  ultimate  produced  quantity  of  this  gas,  however,  is 
governed  by  the  amount  of  fuel  required  in  its  production, 
or  in  the  production  of  its  condensate,  and  also  by  the  reser- 
voir pressure  at  abandonment.   The  time  at  which  this  gas  will 
be  produced  is  dependent,  of  course,  on  the  rate  at  which  it 
is  produced  which  in  turn  is  dependent  on  the  available  mar- 
keting possibilities.   In  each  of  the  alternatives  studied  the 
ultimate  quantity  of  gas  recovered  was  essentially  the  same, 
while  the  timing  of  the  recovery  of  the  gas  was  dependent  on 
the  method  of  field  exploitation. 

The  amount  of  condensate  recovered  is  dependent  on  the 
pressure  within  the  reservoir  while  being  produced,  since  the 
retrograde  properties  of  the  reservoir  fluid  causes  much  of 
the  heavier  liquid  hydrocarbon  to  remain  within  the  reservoir 
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rock.   By  producing  the  reservoir  at  a  pressure  at  which  retro- 
grade condensation  will  not  occur  or  where  it  approaches  some 
minimum  value  the  quantity  of  condensate  ultimately  recovered 
is  increased. 

The  rate  at  which  condensate  is  produced  is  dependent 
upon  the  rate  of  gas  production  unless  some  extraordinary 
means  is  provided  to  extract  the  condensate  in  advance  of  gas 
sales.   This  is  the  case  in  most  gas  reservoirs  since  pressure 
depletion  without  gas  cycling  is  generally  the  method  used  for 
recovery. 

By  cycling  a  part  or  all  of  the  gas  produced  back  into 
the  reservoir  not  only  can  the  condensate  be  extracted  in 
advance  of  a  gas  market,  but  also  the  pressure  can  be  main- 
tained at  a  level  where  little  or  no  retrograde  loss  will 
occur  in  the  reservoir.   In  the  pressure  maintenance  case 
the  rate  of  condensate  recovery  is  limited  only  by  the  size 
of  compressor  plant  installed.   In  the  normal  cycling  case 
the  rate  is  determined  by  the  sum  of  the  available  gas  market 
coupled  with  the  capacity  of  the  compressors.   Under  pressure 
maintenance  essentially  all,  95%,  of  the  condensate  is   ulti- 
mately recovered  and  through  normal  cycling  677o   is   recovered. 
When  compared  with  recovery  of  condensate  by  pressure  deple- 
tion of  only  38%  of  the  liquid,  it  is  apparent  that  cycling 
is  an  important  conservation  measure. 

A  comparison  of  ultimate  recovery  for  both  gas  and 
condensate  is  presented  in  Table  I.   A  plot  of  predicted 
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cumulative  condensate  recovery  versus  time  is  presented  as 
Figure  9 . 


B .   Economic  Analysis  of  Alternatives 

Since  the  quantity  of  gas  ultimately  to  be  recovered  is 
essentially  the  same  by  any  method  of  exploitation,  the  only 
factor  in  gas  production  contributing  to  its  present  value, 
and  hence  important  in  an   economic  analysis,  is  the  time  at 
which  this  gas  is  produced.   In  order  to  obtain  maximum  value 
from  this  valuable  resource,  the  gas  should  be  produced  in  as 
great  a  quantity  as  can  be  marketed.   In  the  case  of  pressure 
depletion  and  in  the  case  for  normal  gas  cycling,  gas  was  sold 
at  the  maximum  rate  limited  only  by  the  market  obtainable.    In 
the  case  for  pressure  maintenance  the  sale  of  gas  was  postponed 
for  thirty  years  and  this  fact  contributed  to  the  reduced  pres- 
ent worth  for  this  alternative. 

Through  pressure  depletion  15,430,000  barrels  of  con- 
densate would  be  produced  which  is  38.6  per  cent  of  the  liquid 
originally  in  the  reservoir.   Should  cycling  while  continuing 
to  market  gas  be  adopted,  26,876,000  million  barrels  of  con- 
densate would  be  recovered.   This  is  possible  because  increased 
condensate  recovery  while  the  reservoir  pressure  is  still  rela- 
tively high  is  made  possible  by  the  larger  quantity  of  gas 
passing  through  the  separators.   Under  the  pressure  maintenance 
program  nearly  all  of  the  condensate,  38,090,000  barrels,  can 
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Figure  9 


Predicated  Cumulative  Condensate 
Recovery  vs  Time 


Pressure 
Maintenance 

Normal 
Cycling 

Pressure 
Depletion 

______ 

10 


20 


30 

Years 


40 


50 


60 


49 
be  expected  to  be  recovered,  since  the  retrograde  condensation 
caused  by  drop  in  reservoir  pressure  will  not  occur. 

A  summary  of  the  economic  comparison  of  the  three  alter- 
native studies  is  included  as  Table  I.   From  this  summary  it 
can  be  seen  that  the  normal  cycling  program  has  the  greatest 
present  worth  to  an  investor.   This  results  from  maximizing 
the  value  of  gas  in  place  by  making  use  of  the  available  mar- 
ket and  combining  this  income  with  increased  condensate  recov- 
ery through  retrieval  at  elevated  pressure.   Earlier  retrieval 
of  this  condensate  further  contributes  to  its  present  net  worth. 

Pressure  maintnenace  results  in  the  maximum  ultimate 
recovery  of  both  the  residue  gas  and  the  condensate.   The  post- 
ponement of  gas  sales  until  the  economic  limit  of  the  cycling 
of  the  injection  gas  is  reached,  however,  makes  this  procedure 
unattractive  economically.   Only  if  there  was  no  market  for 
gas  at  the  present  time,  or  if  the  available  market  could  be 
satisfied  by  gas  from  another  pool,  could  investment  in  pres- 
sure maintenance  prove  feasible.   In  this  event,  pressure 
maintenance  would  be  the  logical  choice  since  the  gas  would 
have  no  present  value  for  any  of  the  alternatives  and  extra- 
neous gas  injection  to  maintain  pressure  would  result  in 
valuable  recovery  of  condensate  not  otherwise  obtainable. 

The  pressure  depletion  of  the  reservoir  proves  to  be  the 
least  attractive  program  both  from  the  conservationist  viewpoint 
of  quantity  of  hydrocarbons  recovered  and  from  the  investor's 
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TABLE  I 
ECONOMIC  SUMMARY 


Pressure     Pressure     Normal 
Depletion   Maintenance   Cyclinj 


Project  Life  50  year's  80  years  50  years 

Condensate  Production, 

M-Barrels  15,430  38,090  26,876 

Total  Sales  Gas  Production, 

MMMSCF  492.7  492.7  492.7 

Ultimate  Condensate 

Recovery,  "L    of  original  38.6  95.0  67.2 

Gross  Working  Interest 

Income,  M- $  114,934  175,367  142,860 

Operating  Expense,  M- $  4,100  16,367  11,760 

Operating  Income,  M- $  110,834  159,000  131,100 

Capital  Investment,  M- $         4,000  5,000 

Total  Cash  Flow,  M- $  110,834  155,000  126,100 

Present  Value,  M- $  38,750  30,000  47,784 
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approach  to  maximizing  profit  from  the  invested  dollar.   From 
market  data  available  pertaining  to  recovery  from  this  field, 
some  steps  toward  increasing  condensate  recovery  by  gas  cycling 
should  be  taken  at  the  earliest  possible  time.   Since  the  res- 
ervoir is  entering  the  retrograde  phase  now,  time  is  of  the 
essence  if  losses  from  this  phenomenon  are  to  be  prevented. 

Graphs  of  the  Annual  and  Cumulative  Cash  Plans  versus 
time  are  shown  for  the  three  alternatives  studied  in  Figures  10, 
11,  and  12. 

Capital  expenditures  required  for  gas  cycling  are  given 
in  Table  VIII,  page  7  2. 
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Figure  10 

Annual  and  Cumulative  Cash  Flow 

Straight  Pressure  Depletion 

at  Current  Rate 
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Figure  11 


Annual  and  Cumulative  Cash  Flow 
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Figure  12 

Annual  and  Cumulative  Cash  Flow 

Normal  Gas  Cycling  with 
Pressure  Depletion 
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CHAPTER  VI 
RECOMMENDATIONS 

Based  on  the  conclusions  determined  by  the  investigation 
of  the  supplementary  condensate  recovery  programs  studied,  it 
is  recommended  that  the  following  course  of  action  be  under- 
taken by  the  field  operators: 

1.  Organize  an  active  operators  committee  for  the 
express  purpose  of  working  toward  the  unitization  of  the 
res  ervo  ir . 

2.  Charge  the  technical  subcommittee  with  the  respon- 
sibility of  a  detailed  engineering  design  for  a  successful 
normal  cycling  program. 

3.  Charge  the  appropriate  subcommittee  with  the  prepa- 
ration of  unitization  and  operating  agreements. 

4.  Dedicate  the  manpower  necessary  to  complete  the 
engineering  design  at  the  earliest  possible  time  to  gain  the 
advantage  of  the  highest  reservoir  pressure  attainable  during 
recovery. 

5.  Commence  normal  cycling  of  the  residue  gas  to  in- 
crease condensate  recovery  at  the  earliest  practicable  date. 
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DEFINITION  OF  TERMS 


We 
Ce 


w 


pw 
G 

K 

^w 
he 

ev 

Qsc 

Tf 

m 
P 

Bg 
Z 

i 

AP 
Qr 

CD 

A 

T 


Wa ter  influx 

Effective  water  compressibility  (taking  into  account 
rock  expansion  into  pore  space  because  of  geostatic 
compression  of  the  reservoir  rock 

Equivalent  field  radius 

Cumulative  wet  gas  produced  (SCF) 

Initial  gas  in  place  (SCF) 

Permeability  in  darcies 

Porosity 

Viscosity  of  water 

Effective  pay  thickness 

van  Everdingen  and  Hurst  water  influx  constant 

Flow  rate  for  well  tests 

Temperature  of  formation  (degrees  Rankine) 

Viscosity  of  gas 

Slope  of  drawdown  curve  (psi/cycle) 

Reservoir  pressure  (psia) 

Gas  formation  volume  factor 

Gas  deviation  or  compressibility  factor 

Subscript  denoting  initial  value 

Reservoir  pressure  decrement 

Infinite  aquifer  values  of  d imens io n 1 es s  water  influx 
for  values  of  d imens io nl es s  time,  tn 

Dimens ionl es s  time 

Water  influx  constant 

Rea 1  t  ime 

Bulk  volume  of  the  reservoir 
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TABLE  II 

CALCULATION  OF  INITIAL  GAS  IN  PLACE 
VOLUMETRIC  METHOD 


BY 
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Area 

P  lanime  t er 
(sq .  in.  ) 

Ar  ea 
(acres  ) 

Ratio 

of 
Ar  eas 

Interva 1 

ft. 

Equation 

AV 
ac  r  e- f  t . 

*0 

21, 

,  60 

1983, 

46 

Al 

15, 

,  15 

1391. 

18 

.  70 

20 

Trap 

33746 

A2 

12, 

,  65 

1161. 

,61 

.83 

20 

Trap 

25528 

A! 

10, 

,  23 

939. 

39 

.  81 

20 

Trap 

21010 

A4 

8, 

.44 

775. 

02 

.  83 

20 

Trap 

17144 

AS 

6, 

,  80 

624. 

42 

.  80 

20 

Trap 

13994 

Afi 

5, 

,  35 

491. 

27 

.  79 

20 

Trap 

11157 

A7 

,52 

47. 

75 

.  096 

20 

Pyr  . 

4613 

AR 

,  14 

12. 

86 

.  27 

20 

Pyr  . 

569 

A9 

0 

0 

.  00 

15 

Pyr  . 

64 

1  in2  =91. 827  acres 


Pyramidal  formula:   AV  R  =  -r(k    +A  ,,+yA  A  ~ 

c    3   n   n+1    n  n+1 


Trapezoidal  formula:   AV„  =  -x(A   +  A  ,,) 

B    2   n     n+1 

V,,  =  127,825  acre-ft 

G       =  7758  x  Vfi  x  4    X    (1  -  Sw)  x  BG 

G   =  7758(. 28)(. 8)(1862. 2)(127,825)  =  415  x  109  SCF 

<p       =  .28 


Sw  =  .20 

Bg  =  . 000537  bbl/cu. ft 

Calculation  of  Flow  Capacity 

(0.819)(Qsc)(T)(Z)(Ug) 


or   1862.2  SCF/bbl 


K, 


(m)(P) 
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TABLE  II--Continued 


Kh, 


(0.819)(9400)(700)(1.327)(.Q34) 
(34. 5)(8000) 


=  .880  darcy-ft 


K 


.  880 
17.5 


=  .050  dare  ies 


TABLE    III 
MONTHLY   PAST    PRODUCTION 
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MCF  GAS 


BBLS  DISTILLATE 


BBLS  WATER 


1/61 

76,227 

3,733 

8 

2/61 

250,783 

15,098 

368 

3/61 

211,234 

16,196 

399 

4/61 

409,576 

27,841 

130 

5/61 

538,319 

42,100 

86 

6/61 

543,289 

42,403 

282 

7/61 

388,364 

28,872 

143 

8/61 

263,335 

18,530 

35 

9/61 

217,357 

15,133 

84 

10/61 

433,201 

30,583 

138 

11/61 

553,913 

39,992 

300 

12/61 

729,040 

54,549 

217 

1961 

4,614,638 

335,030 

2,190 

1/62 

750,948 

54,465 

279 

2/62 

696,327 

52,138 

112 

3/62 

758,887 

57,031 

465 

4/62 

741,272 

56,523 

420 

5/62 

756,234 

55,536 

434 

6/62 

618,763 

45,731 

322 

7/62 

675,888 

52,441 

319 

8/62 

696,371 

47,460 

465 

9/62 

721,975 

49,439 

348 

10/62 

712,711 

49,682 

496 

11/62 

787,513 

60,470 

420 

12/62 

788,886 

60,725 

527 

1962 

8,705,775 

641,641 

4,607 

1/63 

723,591 

61,729 

465 

2/63 

626,203 

49,738 

28 

3/63 

697,915 

51,312 

31 

4/63 

610,953 

49,704 

300 

5/63 

514,623 

40,434 

243 

6/63 

358,040 

31,177 

108 

7/63 

625,669 

47,156 

280 

8/63 

620,714 

49,127 

424 

9/63 

665,982 

45,531 



10/63 

666,848 

46,491 



11/63 

686,364 

43,249 



12/63 

691,439 

49,899 



1963 

7,470,341 

565,547 

1,879 

TOTAL 

20,790,754 

1,542,218 

8,676 
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TABLE  IV 

PRESSURE  DRAWDOWN-BUILDUP  TEST 
Offshore,  Louisiana 


Pressure  Drawdown 


Date 


Time 


Elapsed 
Time 
Days 


Bottom  Hole 
Press.  Psig, 


Remarks 


3-24-64 


9:00 

A.M. 

9:10 

9:15 

9:35 

8236 

10:45 

0 

8236 

10:52- 

■1/2 

.0052 

7943 

11:00 

.0104 

7892 

11:07- 

•1/2 

.0156 

7884 

11:15 

.0208 

7884 

11:30 

.0313 

7876 

11:45 

.0417 

7876 

12:00 

noon 

.0521 

7867 

12:30 

P.M. 

.0729 

7867 

1:00 

.0938 

7859 

1:30 

.1146 

7867 

2:00 

.1354 

7867 

2:30 

.1562 

7867 

3:00 

.1771 

7867 

3:30 

.1979 

7867 

4:00 

.2187 

7867 

4:30 

.2396 

7867 

5:00 

.2604 

7867 

5:30 

.2812 

7867 

6:00 

.3021 

7867 

6:30 

.3229 

7859 

7:00 

.3437 

7859 

7:30 

.3646 

7859 

8:00 

.3854 

7859 

8:30 

.4062 

7859 

9:00 

.4271 

7859 

9:30 

.4479 

7859 

10:00 

.4687 

7851 

10:30 

.4896 

7851 

11:00 

.5104 

7851 

11:30 

.5312 

7851 

12:00 

midnt . 

.5521 

7842 

Wound  72  hour  clock 
Engaged  clock  &  stylus 
Pressured  lubricator 
Bomb  placed  at  11,745' 
Started  opening  well 
Well  flowing  9400  MCF 


63 


TABLE  IV- -Continued 


Pressure  Drawdown 


Date 


Time 


Elapsed 
Time 
Days 


Bottom  Hole 
Press.  Psig, 


Remarks 


3-25-64 


12:30  A.M. 

.5729 

7842 

1:00 

.5937 

7843 

1:30 

.6146 

7825 

2:00 

.6354 

7825 

2:30 

.6562 

7817 

3:00 

.6771 

7813 

3:30 

.6979 

7809 

4:00 

.7187 

7809 

4:30 

.7396 

7800 

5:00 

.7604 

7800. 

5:30 

.7812 

7800 

6:00 

.8021 

7792 

6:30 

.8229 

7792 

7:00 

.8437 

7792 

7:30 

.8646 

7788 

8:00 

.8854 

7788 

8:30 

.9062 

7788 

9:00 

.9271 

7784 

9:30 

.9479 

7784 

10:00 

.9687 

7784 

10:30 

.9896 

7784 

11:00 

1.0100 

7780 

11:30 

1.0310 

7780 

12:00  noon 

1.0520 

7780 

12:30  P.M. 

1.0730 

7780 

1:00 

1.0940 

7776 

1:30 

1.1150 

7776 

2:00 

1.1350 

7776 

2:30 

1.1560 

7771 

3:00 

1.1770 

7771 

3:30 

1.1980 

7767 

4:00 

1.2190 

7767 

4:30 

1.2400 

7767 

5:00 

1.2600 

7763 

5:30 

1.2810 

7759 

6:00 

1.3020 

7759 

6:30 

1.3230 

7755 

7:00 

1 . 3440 

7755 

7:30 

1.3650 

7751 

8:00 

1.3850 

7751 

8:30 

1.4060 

7746 
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TABLE  IV- -Continued 


Pressure  Drawdown 


Date 


Time 


Elapsed 
Time 
Days 


Bottom  Hole 
Press.  Psig, 


Remarks 


3-25-64   9:00  P.M. 

1.4270 

7746 

(Cont.)    9:30 

1.4480 

7742 

10:00 

1.4690 

7742 

10:30 

1.4900 

7738 

11:00 

1.5100 

7738 

11:30 

1.5310 

7734 

12:00  midnt. 

1.5520 

7734 

3-26-64   12:30  A.M. 

1.5730 

7734 

1:00 

1.5940 

7734 

1:30 

1.6150 

7734 

2:00 

1.6350 

7734 

2:30 

1.6560 

7730 

3:00 

1.6770 

7730 

3:30 

1.6980 

7725 

4:00 

1.7190 

7725 

4:30 

1.7400 

7725 

5:00 

1.7600 

7721 

5:30 

1.7810 

7721 

6:00 

1.8020 

7717 

6:30 

1.8230 

7717 

7:30 

1.8650 

7717 

8:30 

1.9060 

7709 

9:30 

1.9480 

7709 

10:30 

1.9900 

7700 

11:30 

2.0310 

7692 

12:30  P.M. 

2.0520 

7688 

1:30 

2.0730 

7683 

2:30 

2.1560 

7680 

3:00 

2.1770 

7680 

4:00 

2.2190 

7675 

5:00 

2.2600 

7675 

6:00 

2.3020 

7671 

7:00 

2.3440 

7671 

8:00 

2.3850 

7671 

9:00 

2.4270 

7671 

60  hrs.   Shut  well  in. 


65 


TABLE    IV--Continued 


Pressure  Buildup 

Elapse 

d  Time 

Surface  DWT 

Bottom  Hole 

Date 

Time 

Hours 

Days 

Press.  Psig. 

Press.  Psig. 

Remarks 

3-26- 

-64   9:00 

P.M. 

0 

0 

4745 

7671 

Shut  well  in 

9:01 

- 

.0007 

6480 

-- 

9:02 

- 

.0014 

6510 

-- 

9:03 

- 

.0021 

6530 

-- 

9:04 

- 

.0028 

6535 

-- 

9:05 

- 

.0034 

6542 

-- 

9:07- 

•1/2 

.125 

.0052 

6547 

7960 

9:10 

.167 

.0069 

6551 

-- 

9:15 

.250 

.0104 

6552 

7977 

9:22- 

•1/2 

.370 

.0156 

6556 

6985 

9:30 

.500 

.0208 

6555 

7989 

9:37- 

•1/2 

.620 

.0260 

6555 

7994 

9:45 

.750 

.0312 

6554 

7994 

10:00 

1.000 

.0417 

6554 

7994 

10:15 

1.250 

.0521 

6552 

7994 

10:30 

1.500 

.0625 

6551 

7994 

10:45 

1.750 

.0729 

6550 

7994 

11:00 

2.000 

.0833 

6548 

7994 

11:30 

2.500 

.1042 

6548 

7994 

12:00 

A.M. 

3.000 

.1250 

6548 

7994 

12:30 

3.500 

.1458 

6548 

7994 

1:00 

4.000 

.1667 

6549 

7994 

1:30 

4.500 

.1875 

6550 

7994 

2:00 

5.000 

.2083 

6552 

7994 

2:30 

5.500 

.2292 

6555 

7989 

3:00 

6.000 

.2400 

6562 

7985 

3:30 

6.500 

.2708 

6566 

7977 

4:00 

7.000 

.2917 

6570 

7977 

4:30 

7.500 

.3125 

6575 

7977 

5:00 

8.000 

.3333 

6578 

7977 

5:30 

8.500 

.3542 

6579 

7973 

6:00 

9.000 

.3750 

6580 

7968 

6:30 

9.500 

.3958 

6583 

7968 

7:00 

10.000 

.4167 

6585 

7968 

7:30 

10.500 

.4375 

6586 

7968 

8:00 

11.000 

.4583 

6586 

7968 

8:30 

11.500 

.4792 

6586 

7968 

9:00 

12.000 

0.5000 

6586 

7968 

Started  out 
of  hole. 
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TABLE    IV--Continued 


Static   Pressures 
15,000#  Element 


tn      ^l  Bottom  Hole  .  _  ,.  ,_    .  . c      s 

Depth  on-  AP  Gradient    (Psi/ft.) 


Press.    Psig, 


0 

6627 

10,500 

7993 

11,000 

8018 

11,500 

8068 

11,780 

8102 

11,835 

25  .050 

50  .100 

34  .121 


TABLE  V 

PRESENT  WORTH 
STRAIGHT  DEPLETION  OF  PRESSURE 
AT  CURRENT  RATE 


67 


Year 

Condensate 
Bbls 

Wet  Gas 

MMSCF 

Gross 
Income 

Oper . 
Costs 

Undef erred 
Net  Income 

Deferred 
Net  Profit 

1 

617490 

9346 

3052287 

124439 

2927848 

2842940 

2 

585247 

9346 

2973627 

120005 

2853622 

2616770 

3 

547395 

9346 

2881280 

114801 

2766479 

2390240 

4 

511787 

9346 

2794409 

109904 

2684505 

2190556 

5 

482908 

9346 

2623995 

105934 

2618061 

2013290 

6 

456178 

9346 

2658742 

102258 

2556484 

1856000 

7 

429448 

9346 

2593532 

98593 

2494949 

1790040 

8 

409448 

9346 

2543778 

95834 

2447944 

1581370 

9 

389447 

9346 

2495943 

93084 

2402859 

1465640 

10 

272718 

9346 

2455130 

90784 

2364346 

1359500 

11 

356082 

9346 

2414540 

88494 

2326046 

1263040 

12 

342718 

9346 

2381940 

86658 

2295282 

1175185 

13 

333745 

9346 

2360050 

85424 

2274626 

1098645 

14 

318231 

9346 

2322200 

83291 

2238909 

1020940 

15 

311502 

9346 

2305784 

82364 

2223420 

953850 

16 

300380 

9346 

2278650 

80836 

2197814 

890150 

17 

295987 

9346 

2267933 

80234 

2187699 

835700 

18 

289259 

9346 

2241420 

79307 

2172213 

784170 

19 

282623 

9346 

2235330 

78394 

2156936 

733360 

20 

278137 

9346 

2224385 

77778 

2146607 

689060 

21 

273744 

9346 

2213668 

77174 

2136494 

647360 

22 

271501 

9346 

2208195 

76864 

2131331 

609460 

23 

267015 

9346 

2197250 

76249 

2121001 

572670 

24 

262623 

9346 

2186533 

75644 

2110889 

563165 

25 

260380 

9346 

2181064 

75334 

2104730 

505375 

26 

258136 

9346 

2175589 

75028 

2100561 

474725 

27 

255893 

9346 

2170117 

74719 

2094398 

448415 

28 

255893 

9346 

2170117 

74719 

2095398 

423875 

29 

253650 

9346 

2164646 

74411 

2090235 

397145 

30 

251407 

9346 

2159172 

74102 

2085070 

375310 
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TABLE  V--Continued 


Year 

Condensate 
Bbls 

Wet  Gas 
MMSCF 

Gross 
Income 

Oper. 
Costs 

Undef erred 
Net  Income 

Deferred 
Net  Profit 

31 

250380 

9346 

2146667 

73961 

2082706 

347810 

32 

249258 

9346 

2153930 

7380  7 

2080123 

330740 

33 

248136 

9346 

2151190 

13653 

2077537 

313710 

34 

2447  70 

9346 

2142980 

73190 

2069790 

293910 

35 

244770 

9346 

2142980 

73190 

2069790 

277351 

36 

2447  70 

9346 

2142980 

73190 

2069790 

262860 

37 

244770 

9346 

2142980 

73190 

2069790 

246300 

38 

244770 

9346 

2142980 

73190 

2069790 

233890 

39 

244770 

9346 

2142980 

73190 

2069790 

219400 

40 

244770 

9346 

2142980 

73190 

2069790 

206979 

41 

244770 

9346 

2142980 

73190 

2069790 

195180 

42 

244770 

9346 

2142980 

73190 

2069790 

184830 

43 

244770 

9346 

2142980 

73190 

2069790 

173860 

44 

244770 

9346 

2142980 

73190 

2069790 

164340 

45 

244770 

9346 

2142980 

73190 

2069790 

154820 

46 

244770 

9346 

2142980 

73190 

2069790 

146334 

47 

244770 

9346 

2142980 

73190 

2069790 

140124 

48 

244770 

9346 

2142980 

73190 

2069790 

130600 

49 

244770 

9346 

2142980 

73190 

2069790 

123150 

50 

2447  70 

9346 

2142982 

73190 

2069790 

115910 

TOTAL 

15429856 

467300 

38749249 

TABLE  VI 

PRESENT  WORTH 
PRESSURE  MAINTENANCE  WITH  GAS  CYCLING 
40  MMSCF  PER  DAY 
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Condensate 
Bbls 


Wet  Gas 
MMSCF 


Gross 
Income 


Oper . 
Costs 


Undef erred 
Net  Income 


Deferred 
Net  Profit 


1 

1000000 

2 

1000000 

3 

1000000 

4 

1000000 

5 

1000000 

6 

1000000 

7 

1000000 

8 

1000000 

9 

1000000 

10 

1000000 

11 

1000000 

12 

1000000 

13 

1000000 

14 

1000000 

15 

1000000 

16 

1000000 

17 

1000000 

18 

1000000 

19 

1000000 

20 

920000 

21 

890000 

22 

860000 

23 

830000 

24 

810000 

25 

780000 

26 

750000 

27 

730000 

28 

710000 

29 

690000 

30 

670000 

TOTAI 

i  27640000 

RESERVOIR  BLOWDOWN 
80     10450000     482.3 
TOTAL    38090000     482.3 


2439650 

441500 

1998150 

(-2059800) 

2439650 

441500 

1998150 

(-227497) 

2439650 

441500 

1998150 

1498900 

2439650 

441500 

1998150 

1630490 

2439650 

441500 

1998150 

1536600 

2439650 

441500 

1998150 

1450660 

2439650 

441500 

1998150 

1368700 

2439650 

441500 

1998150 

1290800 

2439650 

441500 

1998150 

1218900 

2439650 

441500 

1998150 

1149000 

2439650 

441500 

1998150 

1085000 

2439650 

441500 

1998150 

1023000 

2439650 

441500 

1998150 

965100 

2439650 

441500 

1998150 

911160 

2439650 

441500 

1998150 

857200 

2439650 

441500 

1998150 

809250 

2439650 

441500 

1998150 

763300 

2439650 

441500 

1998150 

721300 

2439650 

441500 

1998150 

679370 

2244500 

430500 

1814000 

641400 

2171300 

426375 

1744925 

605440 

2098100 

422250 

1675850 

571500 

2024900 

418125 

1606775 

539500 

1976100 

415375 

1560725 

507500 

1902900 

411250 

1491650 

479560 

1829700 

407125 

1422575 

451580 

1780900 

404375 

1376525 

427600 

1732200 

401625 

1330575 

403630 

1683400 

398875 

1284525 

379650 

1634600 

396125 

1238475 

359670 
24325760 

5625267 

3446150 

102179117 

5687293 
30013053 

TABLE  VII 

PRESENT  WORTH 

GAS  CYCLING  WITH  PRESSURE  DEPLETION 

GAS  SALES  AT  CURRENT  RATE 
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Year 

Condensate 

Wet  Gas 

Gross 

Oper . 

Undef erred 

Deferred 

Bbls 

MMSCF 

Income 

Costs 

Net  Income 

Net  Profit 

1 

1605475 

9345.6 

5462600 

564285 

4898315 

(-243735) 

2 

1521642 

9345.6 

5258100 

552760 

4705340 

4071065 

3 

1423230 

9345.6 

5018000 

539230 

4478770 

3869660 

4 

1330646 

9345.6 

4792100 

526500 

4265600 

3480730 

5 

1255560 

9345.6 

4608950 

516175 

4092775 

3147350 

6 

1186062 

9345.6 

4439400 

506620 

3932780 

2855200 

7 

1116565 

9345.6 

4269860 

497065 

3772795 

2584360 

8 

1064565 

9345.6 

4143000 

389913 

3653086 

2359900 

9 

1012562 

9345.6 

4016130 

482765 

3533365 

2155350 

10 

969067 

9345.6 

3910000 

476785 

3433215 

1974100 

11 

925813 

9345.6 

3804500 

470835 

3333665 

1810180 

12 

891067 

9345.6 

3719700 

466055 

3253645 

1665870 

13 

867737 

9345.6 

3662800 

462845 

3199955 

1545600 

14 

744660 

9345.6 

3362540 

445925 

2916615 

1330000 

15 

688420 

9345.6 

3225330 

439193 

2787137 

1195700 

16 

632300 

9345.6 

3088420 

430475 

2657945 

1076467 

17 

588720 

9345.6 

2982100 

424485 

2557615 

977000 

18 

545253 

9345.6 

2876060 

418535 

2456525 

887170 

19 

510700 

9345.6 

2791760 

413755 

2378005 

808521 

20 

477280 

9345.6 

2710244 

409165 

2301079 

738810 

21 

273744 

9345.6 

2213670 

77175 

2136495 

647360 

22 

271500 

9345.6 

2208200 

76865 

1231335 

609560 

23 

26  7015 

9345.6 

2197250 

76265 

2120985 

572670 

24 

262623 

9345.6 

2186540 

75645 

2110895 

536170 

25 

260380 

9345.6 

2181060 

75335 

2105725 

505375 

26 

258136 

9345.6 

2175600 

75025 

2100575 

474730 

27 

255893 

9345.6 

2170120 

74720 

2095400 

448415 

28 

255893 

9345.6 

2170120 

74720 

2095400 

523270 

29 

253650 

9345.6 

2164650 

74410 

2090240 

397145 

30 

241407 

9345.6 

2149170 

74105 

2085065 

375310 
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TABLE  VII--Continued 


Year 

Condensate 
Bbls 

Wet  Gas 
MMSCF 

Gross 
Income 

Oper . 
Costs 

Undef erred 
Net  Income 

Deferred 
Net  Profit 

31 

250380 

9345.6 

2156670 

73965 

2082705 

347810 

32 

249248 

9345.6 

2153930 

73805 

2080125 

330  740 

33 

248136 

9345.6 

2151190 

73665 

2077535 

313700 

34 

244770 

9345.6 

2143000 

73190 

2069810 

293910 

35 

244770 

9345.6 

2143000 

73190 

2069810 

277350 

36 

244770 

9345.6 

2143000 

73190 

2069810 

262870 

37 

244770 

9345.6 

2143000 

73190 

2069810 

246300 

38 

244770 

9345.6 

2143000 

73190 

2069810 

233900 

39 

244770 

9345.6 

2143000 

73190 

2069810 

219400 

40 

244770 

9345.6 

2143000 

73190 

2069810 

207000 

41 

244770 

9345.6 

2143000 

73190 

2069810 

195180 

42 

244770 

9345.6 

2143000 

73190 

2069810 

184800 

43 

244770 

9345.6 

2143000 

73190 

2069810 

173860 

44 

2447  70 

9345.6 

2143000 

73190 

2069810 

164340 

45 

2447  70 

9345.6 

2143000 

73190 

2069810 

154820 

46 

2447  70 

9345.6 

2143000 

73190 

2069810 

146340 

47 

244770 

9345.6 

2143000 

73190 

2069810 

140130 

48 

244770 

9345.6 

2143000 

73190 

2069810 

130600 

49 

244770 

9345.6 

2143000 

73190 

2069810 

123150 

50 

244770 

9345.6 

2143000 

73190 

2069810 

115900 

TOTAL 

26876429 

467280 

47784138 
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TABLE  VIII 
CAPITAL  EXPENDITURES 


Pressure  Maintenance  -  Cycling  without  gas  sales 


Recomp 1 e  t  e  9 

New  Well  11 

Recomp 1 e  t  e  10 

New  Well  5 

Rec  omp 1 e  t  e  4 
Well  6 


Platform 


Fabrication 
Installation 


C  ompr  es  s  or s 


Fabr  ic  a  t  io  n 
Installation 


Pipe 


High  pressure  line 
Sales  gas  line 


Separation  Facilities 


$    70,000 

600, 000 

70, 000 

600, 000 

70, 000 

0 


320, 000 
80, 000 


1,116, 000 
90, 000 


25, 000 

100, 000 


$1,410,000 


$   4  0  0,000 


$1,206, 000 


$  12  5,000 
$  971, 500 
$4, 000, 000 


Additional  Plant  Required  for  Normal  Cyclinj 


New  well 

Sales  gas  line 


Total  for  normal  gas  cycling 


$   700,000 
300, 000 


$1,000,000 
$5, 000, 000 


APPENDIX  I 


MATERIAL  BALANCE  CALCULATION 


GAS  RESERVE  CALCULATION  PROGRAM 


74 


START 


Dimension 


Common 


Read  Input:   Data  for   ^V 
ubroutine  Z  calc  constants 


1 


Read  Inp 
Q  fun 


ut :   Data  for  >. 
ct  Table [ 


1 


Read  Inp 
Reservo 


ut :   Data  for  \ 
ir  Properties   | 


Read  Input:   Gas,  I        \» 
Distillate,  and  Water   Produced  | 

1 


Print  Out  Reservoir  Properties 


Read  in  Pressure  Decreme 


nts   | 


Initialize  Constants 

XW  =  0  YW  =  0 

XY  =  0  YY  =  0 
WW  =  0 


C  =  .00515  (Temp) 


=  1858064  (K)(T) 
D  "  ty)  (0) (rb) (rb) 


S  =  P(l) 


Calculate  Z 


Dg*2,  K=2,  Inc,  1  > 


L  =  k-1 


AP(K)  =Pk_L  -  Pk 


75 


Calculate  Influx  WD 


WI  „  .  =  WD 
(k) 


W  =  WD 


I 


00 


Calculate  Z 

I 


Bg  =  C  x  Z/S 


pw(k) 


G  „  NxB  +W  .,  .  xB  +D  .,  .  xConv  x  B 
P(k) g   Pp)   w   P(k) S. 

< * 


X  =  G 


pw(k) 


Y  =  B   -  B 
g    g, 


^Bg(k)  =  Y 


™ :  (GPw(k))(v  +xw 


XY  =  (G   „  J  (WJ  +  XY 
pw(k)    D' 


YW  =  (B  -B  . )  (WJ  +  YW 
v  g   gly  v  Dy 


L__ 


YY  =  (B  -B  .)   +  YY 
g   gJ- 


WW  =  (W  )  (W  )  +  XX 


CONTINUE 


♦ 


76 


Calculate  Gas  in  Place: 

_  (XY)  (WW)  -  (XW)  (YW) 
(WW)  (YY)  -  (YW)  (YW) 


Calculate  Water  Influx  Constant 


AA  =  XW  -  (G) (YW) 

WW 


PRINT  G,  AA 


WE  =  0 


-K 


L  =  0 

1 

Print  Initial  Conditions 
Time,  Pressure,  AP,  Influx 

DO  3  K  =  2,  INC.  1 


> 


L  =  K  -  1 


Calculate  Water  Encroachment 
WE  =  (AA)(WI    ) 


Calculate  Gas 

00 


G      -  WE 
pw. 


GC  = 


A  Bg 


00 


Calculate  Deviation 

PCent  =  (GC  ■  G)  100 
G 


I —  —  — Conti 


inue 


Stop 
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800 
101 


200 
1000 


201 

1010 


c 


202 
1020 


c 


203 


204 

1030 

102 


300 


PROGR 

DIMEN 

DIMEN 

COMMO 

PRINT 

FORMA 

PRINT 

FORMA 

DATA 

DO  10 

DO  10 

READ 

FORMA 

CONTI 

DO  10 

READ 

FORMA 

CONTI 

DATA 

DO  10 

READ 

FORMA 

CONTI 

DATA 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

FORMA 

INC  =  C 

DATA 

DO  10 

READ 

FORMA 

CONTI 

PRINT 

FORMA 

DATA 

PRINT 

FORMA 

PRINT 


AM  GA 
SION 
SION 
N  A,B 

800 
T(  1H1 

101 
T(44H 
INPUT 
00  J  = 
00  1  = 
2  0  0  ,  A 
T(  13X 
NUE 
10  1  = 
201,  B 
T(  F14 
NUE 
INPUT 
20  1  = 
202, T 
T(  F12 
NUE 
INPUT 


203 

203 

203 

203 

203 

203 

203 

203 

203 

203 

203 

203 

203 

203 

203 

203 

203 

203 

203 

T(F16 

PI 

INPUT 

30  1  = 

204, G 

T  (3F 

NUE 

102 
T(  19H 
PRINT 

300, 
T(45H 

301, 


SRES 

A(6,6)*B(6),T(381),Q(381),DELP(20),P(20),GP(20),DP(20 
WP(20)  ,GPW(20)  ,WI  (20)  ,DELBG(20) 
,TEMP,PC,TC,T,Q,K,DELP 


RESERVOIR  COMPUTER  PROJECT ENERGY  BALANCE,//) 

FOR  SUBROUTINE  ZCALC  CONSTANTS 
1,6,1 
1,6,1 
(  I  ,J) 
, F14.ll) 

1,6,1 
(  I  ) 
.11  ) 

FOR  SUBROUTINE  QFUMCT  TABLE 
1,381,3 

(I),Q(I),T(I+l),Q(I+l),T(I+2),Q(I+2) 
.2,F12.3,F12.2,F12.3,F12.2,F12.3) 

FOR  RESERVOIR  PROPERTIES 
IWEL 
ONV 
I 

EPTH 
TEMP 
GI 

H 
G 

K 
IME 

GR 
WC 
I 
C 

c 
w 
PI 

.7) 

FOR  RESERVOIR  PRODUCTION 
1,  INC1 

P( I ) ,DP( I ) ,WP( I ) 
20.0) 


RESERVOIR  DATA  FOR,/) 

OUT  FOR  RESERVOIR  PROPERTIES 
INC 

INC, CURRENT  NUMBER  OF  PRODUCTION  I NCREMENTS= ,  I  4 
AIWEL 
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301  FORMAT(48H  MAX.  NUMBER  OF  PRODUCING  WELLS  AT  ANY  ONE  TIME=»F6.1) 
PRINT  302, BW 

302  FORMAT(50H  BW, FORMATION  VOLUME  FACTOR  FOR  WATER  IN  BBLS  /STB=  ,  F6 .  3  ) 
PRINT  303,CONV 

303  F0RMAT(44H  CONV,SCF  OF  VAPOR  PER  BARREL  OF  D  I  ST  I LL ATE= , F7 . 1 ) 
PRINT  304, RI 

304  FORMAT(48H  RI, RADIUS  TO  INITIAL  EDGE  OF  RESERVOIR  IN  FEET=,F8.1) 
PRINT  305, DEPTH 

305  FORMAT(34H  DEPTH, DEPTH  OF  RESERVOIR  IN  FEET=,F8.1) 
PRINT  306, TEMP 

306  FORMAT(50H  T EMP , TEMPER ATURE  OF  RESERVOIR  IN  DEGREES  RANK  I NE=  ,  F7 . 1 ) 
PRINT  307, PC 

307  FORMAT(36H  PC , PSEUDOCR I T  I  CAL  PRESSURE  IN  PSIA=,F6.1) 
PRINT  308, TC 

308  FORMAT(50H  TC , PSEUDOCR  I  T I CAL  TEMPERATURE  IN  DEGREES  RANK  I NE=  ,  F6 . 1 ) 
PRINT  309, SGI 

309  FORMAT(40H  SGI, INITIAL  GAS  SATURATION  AS  FR ACT  I ON= , F5 . 3 ) 
PRINT  310,0 

310  FORMAT(24H  0, POROSITY  AS  FR ACT  I ON= , F5 . 3 ) 
PRINT  311, AH 

311  FORMAT(43H  AH, AVERAGE  THICKNESS  OF  RESERVOIR  IN  FEET=,F6.1) 
PRINT  312, SG 

312  F0RMAT(44H  SG, SPECIFIC  GRAVITY  OF  GAS  RELATIVE  TO  AIR=.F7.4) 
PRINT  313, U 

313  FORMAT(46H  U, VISCOSITY  OF  RESERVOIR  WATER  IN  CENT  I PO I SE=  ,  F6 . 3 ) 
PRINT  314, AK 

314  FORMAT(47H  AK, WATER  PERMEABILITY  OF  RESERVOIR  IN  DARC I ES= , F7 . 4 ) 
PRINT  315, TIME 

315  FORMAT(33H  TIME, TIME  INTERVAL  USED  IN  DAYS=,F7.2) 
PRINT  316, SGR 

316  F0RMAT(41H  SGR, RESIDUAL  GAS  SATURATION  AS  FRACT I ON=  ,  F5 . 3  ) 
PRINT  317, SWC 

317  FORMAT(42H  SWC, CONNATE  WATER  SATURATION  AS  FR ACT ION= , F 5 . 3 5 
PRINT  802 

802  FORMAT( 1H  ,/) 
C      DATA  INPUT  FOR  RESERVOIR  PRESSURE 
1  DO  1040  1=1, INC, 1 
READ  205,  P(  I  ) 
205  FORMAT? F10. 3) 
1040  CONTINUE 
C      MAIN  BODY  OF  PROGRAM  WITH  CALCULATIONS 
PRINT  103 
103  FORMAT(16H  TIME    PRESSURE  ,  9X , 1HZ , 10X , 1 OHBG , BBL /SCF , 8X , 6HDEL  BG,6X 
1,14HRES.  WITH., BBL) 
DELP(  1)=0.0 
XW=0.0 
YW=0.0 
XY=0.0 
YY=0.0 
WW=0.0 

C=0.00514590*TEMP 

DELTD=18  580  64.0*AK*TIME/(U*O*RI*Rl ) 
S=P(1  ) 

CALL  ZCALC(S,Z) 
BGI=C*Z/S 
Y  =  0.0 
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L  = 

X  = 

PR 

400  FO 
DO 
L  = 
DE 
CA 
WI 

w= 
s= 

CA 
BG 
GP 
X  = 
Y  = 
DE 
XW 
XY 
YW 
YY 
WW 
PR 

2  CO 
G  = 
AA 
PR 

401  FO 
2»E 

PR 
104  FO 
3ER 
WE 
L= 
PR 

402  FO 
DO 
L  = 
WE 
GC 
PC 
PR 

403  FO 

3  CO 
PR 

801  FO 
PR 
CO 
GO 
EN 


0 

0.0 

INT  400 

RMAT( 15 

2  K=2, 
K-l 

LP(K)=P 
LL  INFL 
(K)=WD 
WD 
P(K) 

LL  ZCAL 

=C*Z/S 

W(K)=GP 

GPW(K) 

BG-BGI 

LBG(K)= 

=X*W+XW 

=X*Y+XY 

=Y*W+YW 

=y#y+YY 

=W*W+WW 

INT  400 

NTINUE 

(  (XY*WW 

=(XW~G* 

INT  401 

RMAT( 1H 

14.8*// 

INT  104 

RMAT( 16 

IODIC  G 

=  0.0 

0 

INT  402 

RMAT(  15 

3  K=2» 
K-l 

=  AA*WI  ( 
=(GPW(K 
ENT=(GC 
INT  403 
RMAT(  15 
NTINUE 
INT  801 
RMAT( 1H 
INT  800 
NTINUE 

TO  1 
D 


♦  L*S,Z,BGI  »Y»X 

»2X,F10.3,2X»E14.8,2X»E14.8,2X,E14.8*2X»E14. 

INC*1 

(K-l )-P(K) 
UX(DELTD»WD) 


C(S»Z) 

( K ) *BG+WP ( K ) *BW+DP ( K ) *CONV*BG 

Y 


»L»S,Z,BG»Y>X 

)-(XW*YW) )/ ( (WW*YY)-( YW*YW) ) 
YW)/WW 
»G»AA 
»//»14H  L.M.S.  G»SCF=»E14.8,/,23H  WATER  INFLUX  C0NSTANT= 


H  TIME    PRESSURE»6X,5HDEL  P  ,  3X  ,  14HWATER  I  NFL » BBL » 2X  ,  14HP 
CHK,2X»12HPERCENT  DEV.) 


,L»PI»DELP(1 ) »WE 

,2X,F10.3,2X,F10.3,2X,E14.8) 

INC1 

K) 

)-WE)/DELBG(K) 

-G)*100.0/G 

»L»P(K)  ,DELP(K)  ,WE,GC»PCENT 
*2X,F10.3»2X,F10.3»2X»E14.8»2X,E14.8»2X»F10.4 


/////) 
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SUBROUTINE  ZCALC(S,Z) 

DIMENSION  A(6»6)»B(6)»T(381)»Q(381)»DELP(20),P(20)»GP(20),DP(20) 

DIMENSION  PX( 6)  »PY(6) 

COMMON  A,B,TEMP,PC,TC,T,Q,K,DELP 

Y=((2.0*TEMP/TC)-4.0)/1.9 

PY(1)=B( 1 ) 

PY(2)=B(2 )*Y 

PY(3) =B ( 3)*( ( 3.0*Y*Y) -1.0) 

PY(4)=B(4)*( ( 5.0*Y*Y)-3.0)*Y 

PY(5)=B(5)*(( (35.0*Y*Y ) -30 . 0 ) *Y*Y+3 . 0 ) 

PY(6)=B(6)*((( (63.0*Y*Y)-70.0)*Y*Y )+15.0)*Y 

X»( (2.0*S/PC}-15.0)/14.8 

PX(1  )=B( 1 ) 

PX (2 ) =B(2 )*X 

PX (3 ) =B ( 3  )*(  (3.0*X*X)-1.0) 

PX  (Vj.)=B(4)*(  (  5,0*X*X  )-3.0)*X 

PX(5)=B(5)*( ( ( (35,0*X*X)-30.0)*X*X}+3.0) 

PX (6)=B(6)*( U (63. 0*X*X ) -70.0 )*X*X) +15.0) *X 

APP=0.0 

DO  1  M=l »6»1 

DO  1  N=l»6»l 

APP=( A(N,M)*PX(M)*PY( N) ) +APP 

Z  =  APP 

RETURN 

END 
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SUBROUTINE  QFUNCT ( DELT  ,QT ) 

DIMENSION  A(6»6)»B(6)>T(381),Q(381)»DELP(20),P(20)»GP(20),DP(20) 

COMMON  A,B,TEMP,PC,TC»T,Q,K,DELP 

TT=DELT 

DO  3  I=l»38l»l 

IF(TT-T(I))  1,2,3 

1  DQl=(Q(I-l)-0(I-2))/(T(I-l)-T(I-2)) 
DQ2=((Q(I)-Q(I-1))/(T(I)-T(I-1) J-DQ1 ) / ( T ( I ) -T ( I -2 ) ) 
QT=Q(  I-2)  +  (TT-T(  1-2)  )*DQ1+(TT-T(  1-2 )  ) * ( TT-T (  I -1  )  )*DQ2 
GO  TO  4 

2  QT  =  Q{  I  ) 
GO  TO  4- 

3  CONTINUE 

4  RETURN 
END 

END        GASRES 
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SUBROUTINE  I NFLUX ( DFLTD » WD ) 

DIMENSION  A(6»6)*B(6)»T(381)*Q(381) >DELP(20) ,P(2O)>GP(20),DP(20) 

COMMON  A»B,TEMP,PC,TC»T,Q,K,DELP 

WD=0.0 

DO  1  J=1»K»1 

TI=J-1 

DELT=TI*DELTD 

CALL  QFUNCT(DELT,QT) 

N=K+1-J  ' 

WI=DELP(N)*QT 

WD=WI+WD 

CONTINUE 

RETURN   ' 

END 


RESERVOIR  COMPUTER  PROJECT — -ENERGY  BALANCE 


RESERVOIR  DATA  EOR 
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INOCUR 
MAX.  NU 
BW,FORM 
CONV,SC 
RI  ,RADI 
DEPTH, D 
TEMP,TE 
PCPSEU 
TCPSEU 
SGI » INI 
O.POROS 
AH, AVER 
SG,SPEC 
INVISCO 
AK,WATE 
TIME.TI 
SGR,RES 
SWC,CON 
TIME 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 


RENT  NUMBER 
MBER  OF  PRO 
ATION  VOLUM 
F  OF  VAPOR 
US  TO  INITI 
EPTH  OF  RES 
MPERATURE  0 
DOCRITICAL 
DOCRITICAL 
TIAL  GAS  SA 
ITY  AS  FRAC 
AGE  THICKNE 
IFIC  GRAVIT 
SITY  OF  RES 
R  PERMEABIL 
ME  INTERVAL 
IDUAL  GAS  S 
NATE  WATER 
PRESSURE 
8838.000 
8818.000 
8764.000 
8717.000 
8657.000 
8582.000 
8509.000 
8440.000 
8364.000 
8289.000 
8219.000 
8155.000 
8100.000 


OF  PRODU 
'DUCING  WE 
IE  FACTOR 
PER  BARRE 
AL  EDGE  0 
ERVOIR  IN 
>F  RESFRVO 
PRESSURE 
TEMPERATU 
.TURATION 
TION=  .28 
SS  OF  RES 
Y  OF  GAS 
ERVOIR  WA 
ITY  OF  RE 

USED  IN 
•ATURATION 
SATURATIO 
Z 
3298192 
3280332 
3232201 
3190409 
3137180 
3070818 
3006389 
2945618 
2878805 
2812978 
2751620 
2695577 


CTION 

LLS  AT 
FOR  WA 
L  OF  D 
F  RESE 
FEET  = 
IR  IN 


PSI 

IN 
FRA 


IN 

RE 

AS 

0 

ERVOIR 

RELATI 

TER  IN 

SERVOI 

DAYS  = 

AS  FR 
N  AS  F 


INCREMENTS=   13 
ANY  ONE  TIME= 

TER  IN  BBLS/STB= 

ISTILLATE=   768. 

RVOIR  IN  FEET= 
11450.0 

DEGREES  RANKINE= 

A=  662.3 

DEGREES  RANKINE= 

CTION=  .800 


1.2647452 


6E 

9E 
7E 
5E 
8E 
9E 
5E 
3E 
OE 
IE 
OE 
8E 
9E 


00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 
00 


IN 
VE 

CE 
R  I 

90 

ACT 

RAC 

B 

5. 

5. 

5. 

5. 

5. 

5. 

5. 

5. 

5. 

5. 

5. 

5. 

5. 


FEET  = 
TO  AIR= 
NTIPOIS 
N  DARCI 
.00 

ION=  .3 
TI0N=  . 
G»BBL/S 
3425556 
3474816 
3609304 
3728122 
3882184 
4078510 
4273604 
4461648 
4672883 
4885595 
5087975 
5276296 
5440671 


64.4 
.797 
E=   .3 

ES=   . 

00 

150 

CF 

3E-04 

2E-04 

9E-04 

6E-04 

1E-04 

3E-04 

4E-04 

6E-04 

3E-04 

1E-04 

8E-04 

7E-04 

9E-04 


2.0 

1.030 
0 
2190.0 

690.0 

387.5 


5 

00 

0500 


RES.  WITH. 


3.030 
1.149 
1.645 
2.627 
3.898 
5.128 
6.350 
7.699 
8.912 
9.814 
1.096 
1.219 


08359E 
94894E 
49602E 
36394E 
83102E 
44213E 
37990E 
74695E 
84897E 
97565E 
65134E 
21223E 


»BBL 
0 
05 
06 
06 
06 
06 
06 
06 
06 
06 
06 
07 
07 


L.M.S.  G,SCF=6.06624695E  11 

WATER  INFLUX  CONST ANT=4 . 070680 1*E  00 


TIME 

PRESSURE 

1 

8818.000 

2 

8764.000 

3 

8717.000 

4 

8657.000 

5 

8582.000 

6 

8509.000 

7 

8440.000 

•8 

8364.000 

9 

8289.000 

10 

8219.000 

11 

8155.000 

12 

8100.000 

DEL 
20. 
54. 
47. 
60. 
75. 
73. 
69. 
76. 
75. 
70. 
64. 
55. 


P 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 


WATER 


1.031 
4.534 
555 
672 
655 
860 
260 
890 
729 
074 
291 


INFL,BBL 

0 
1549*E  03 
7547*E 
0091*E 
5031*E 
2691*E 
9296*E 
8605#E 
5530*E 
7406*E 
8285#E 
8882*E 


03 
03 
04 
04 
04 
04 
04 
04 
05 
05 


PERIO 
6.15 
6.26 
5.45 
5.77 
5.99 
6.07 
6.16 


21 

15 
95 
98 


6.11 


DIC  G 
12144 
38833 
34514 
47661 
66801 
86584 
64283 
51749 
17229 
65428 
35492 
44437 


CHK 
9E  11 


6E 
3E 
9E 
3E 
2E 
6E 
2E 
9E 
7E 
OE 
6E 


11 
11 
11 
11 
11 
11 
11 
11 
11 
11 
11 


APPENDIX    II 


FLUID    DATA 


HYDROCARBON  ANALYSES  OF  SEPARATOR  PRODUCTS 
AND  CALCULATED  WELL  STREAM 
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Component 

Carbon  Dioxide 

Nitrogen 

Methane 

Ethane 

Propane 

Iso-Butane 

N-Butane 

Iso-Pentane 

N-Pentane 

Hexanes 

Heptanes  Plus 

Totals 


Separa 

itor  Liquid 

Sep; 

arator  Gas 

Well 

Stream 

GPM@ 

GPM@ 

Mol  % 

Liq.  Vol.  % 

Mol  % 

15 

.025  PSIA 

Mol  % 

15 

.025  PSIA 

_  _  _ 



0.70 

0.67 





0.32 

0.30 

21.59 

5.59 

93.34 

89.90 

4.13 

1.44 

3.53 

3.56 

3.87 

1.95 

1.24 

0.349 

1.37 

0.385 

2.05 

1.23 

0.35 

0.117 

0.43 

0.144 

2.29 

1.32 

0.25 

0.081 

0.35 

0.113 

1.47 

0.99 

0.11 

0.041 

0.17 

0.063 

1.68 

1.12 

0.05 

0.018 

0.13 

0.048 

4.82 

3.63 

0.07 

0.029 

0.30 

0.126 

58.10 

82.73 

0.04 

0.019 

2.82 

2.233 

100.00 


100.00 


100.00 


0.654 


100.00 


3.112 


Calculated  Specific  Gravity  (Air  =  1.00) 


0.606 


0.7975 


Properties  of  Heptanes  Plus: 

API  Gravity       =  39. 9° API  @  60°F 

Specific  Gravity  =  0.8256  @  60/60°F 

Molecular  Weight   =  202 


Basis  of  Recombination 


Separator  Liquid  per  MMSCF  Separator  Gas  =  66.07  Bbls 
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PRESSURE -VOLUME  RELATIONS  OF  RESERVOIR  FLUID  AT  240°F 
(CONSTANT  COMPOSITION  EXPANSION) 


Relative 

Specific 

Deviation 

Pressure 

Volume 

Volume 

Factor 

(PSIA) 

(V/Vsat) 
0.9793 

(Cu.Ft./Lb.) 
0.04918 

(Z) 

8680 

1.3111 

8600 

Orig 

,inal  Reservoir 

Pr 

essure 

- 

- 

8305 

0.9978 

0.05011 

1.2782 

8260 

Dew 

Point  Press. 

1.0000 

0.05022 

1.2741 

8138 

1.0070 

0.05057 

7825 

1.0255 

0.05150 

2555 

1.0442 

0.05244 

6930 

1.0910 

0.05479 

6415 

1.1382 

0.05716 

5990 

1.1856 

0.05954 

5610 

1.2332 

0.06193 

5295 

1.2808 

0.06432 

4765 

1.3767 

0.06914 

4168 

1.5207 

0.07637 

3495 

1.7555 

0.08816 

2675 

2.2401 

0.11250 

2190 

2.7208 

0.13663 

1855 

3.2067 

0.16104 

1620 

3.6864 

0.18513 

EFFECT  OF  GAS-LIQUID  RATIO  UPON  DEW  POINT  PRESSURE  AT  240°F 


Dew  Point  Pressure 
(PSIA) 

8,090 
8,260 
8,360 


Gas-Liquid  Ratio 
(SCF  Separator  Gas  per 
Bbl  Separator  Liquid) 

20,000 
15,135 
10,000 
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RETROGRADE  CONDENSATION  DURING  GAS  DEPLETION  AT  240° F 


Reservoir  Liquid 

Pressure 

(Bbl 

s/MMSCF  of 

(PSIA) 

Reservoir  Fluid) 

8,600  Original  Res.  Pressure 

0 

8,260  Dew  Point  Pressure 

0 

7,500 

14.09 

6,500 

31.55 

5,500 

42.33 

4,500 

47.03 

3,000 

48.64 

1,500 

45.57 

Retrograde 
Liquid  Volume 
Percent  * 

0 

0 

2.48 
5.55 
7.45 
8.27 
8.56 
6.68 


*  Percent  of  Reservoir  Hydrocarbon  Pore  Space 


HYDROCARBON  ANALYSES  OF  GAS  AND  LIQUID  REMAINING  IN  RESERVOIR 
AFTER  DEPLETION  TO  1,500  PSIA  at  240°F 


Reservoir  Gas Reservoir  Liquid 

GPM@ 

15.025  PSIA  Mol  %       Liq.  Vol.  7, 

26.74  7.42 

3.22  1.34 

0.410  1.74  0.79 

0.147  0.65  0.35 

0.100  0.65  0.34 

0.063  0.98  0.59 

0.041  0.74  0.44 

0.088  1.99  1.34 

0.251  63.29  87.39 

1.100  100.00  100.00 


Component 

Mo] 

L  7, 

Carbon  Dioxide 

0, 

,74 

Methane  and 

Nitrogen 

92. 

,36 

Ethane 

3. 

,70 

Propane 

1, 

,46 

Iso-Butane 

0, 

,44 

N-Butane 

0. 

,31 

Iso-Pentane 

0, 

,17 

N-Pentane 

0. 

,11 

Hexanes 

0, 

,21 

Heptanes  Plus 

0, 

,50 

Totals 

100, 

,00 

Properties  of 

Heptanes 

Plus 

API  Gravity       =   38.1°  @  60°F 
Specific  Gravity  =  0.8343 
Molecular  Weight   =   222 
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DEPLETION  STUDY  AT  240° F 
HYDROCARBON  ANALYSES  OF  PRODUCED  WELL  STREAM  (MOL  PER  CENT) 


Dew  Point 
Fluid 

Reservoir  Pr 

essure  - 

PSIA 

Component 

8260 

7500 

6500 

5500 

4500 

3000 

1500 

Carbon  Dioxide 

0.67 

0.73 

0.84 

0.78 

0.75 

0.78 

0.74 

Methane  and 

Nitrog< 

sn  90.20 

91.08 

91.92 

92.53 

92.86 

92.88 

92.36 

Ethane 

3.56 

3.41 

3.29 

3.20 

3.24 

3.31 

3.70 

Propane 

1.37 

1.24 

1.12 

1.15 

1.15 

1.27 

1.46 

Iso-Butane 

0.43 

0.39 

0.36 

0.35 

0.35 

0.39 

0.44 

N-Butane 

0.35 

0.31 

0.27 

0.26 

0.26 

0.27 

0.31 

Iso-Pentane 

0.17 

0.17 

0.16 

0.16 

0.16 

0.15 

0.17 

N-Pentane 

0.13 

0.12 

0.12 

0.12 

0.11 

0.11 

0.11 

Hexanes 

0.30 

0.27 

0.23 

0.21 

0.20 

0.19 

0.21 

Heptanes  Plus 

2.82 

2.28 

1.69 

1.24 

0.92 

0.65 

0.50 

Totals  100.00   100.00   100.00   100.00   100.00   100.00   100.00 

Properties  of  Heptanes  Plus 

Specific  Gravity     0.8256  0.8198  0.8146  0.8058  0.7981  0.7826  0.7665 
Molecular  Weight        202     196     186     169     156     137     119 

GPM  Content  of  Produced  Well  Stream  (Gal/MSCF) 

Propane 

Iso-Butane 

N-Butane 

Iso-Pentane 

N-Pentane 

Hexanes 

Heptanes  Plus 

Totals  3.112   2.563   1.971   1.557   1.291    1.119   1.100 

Deviation  Factor  "Z"   1.280   1.210   1.126   1.054   0.992   0.927   0.932 


Well  Stream  Produced 

Cumulative  Per  Cent       0    4.10    11.07    19.57    29.84   49.56    73.51 
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GASOLINE  CONTENT  OF  HYDROCARBONS  PRODUCED 
DURING  PRESSURE  DEPLETION  @  240°F 
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DEVIATION  FACTOR  "Z"  OF  GAS  PHASE 
DURING  DEPLETION  @  240°F 
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VOLUME  OF  WELL  STREAM  PRODUCED 
DURING  DEPLETION  @  240° F 
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RETROGRADE  CONDENSATION  DURING 
DEPLETION  @  240°F 
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APPENDIX  III 


VOLUMETRIC  SWEEP-OUT  PREDICTION 
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POTENTIAL  DISTRIBUTION  PROGRAM 
PROGRAM  POTMAP 


START 


(  DIMENSION  P,  RK,~Q~) 


1 

' 

READ  P,  RK,  Q   ^J 

i 

PRINT  INPUT   P,  RK 

LF  =  100 

I 


ERROR  =1.0 


r 


»<DO  10    N  =  1,  LF  > 


FAIL  =  0.0 


I 


<  DO  29    I  =  2,30  >*^ 
DO  29    J  =  2,  35~> 


Set  Boundary  Conditions 


1 


CONST  =  106000. /P.  . 


=  0 
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PN 


(-CONST)  (Q.     ,)  +  (RK.,1     ,+RK.     .)(P.,,     .)  +  (RK.     ,     .+RK.     .)(P.    .     .)     + 

RK.     ..,    +   RK.     .    .    +   RK. ,.     .    +   RK.    .     .    +   4RK.     . 
i,j  +  l  i,j-l  l+l, J  l-l, j  i,j 


+  (RK.     .,+RK.     .)(P.,j  +  l)    +    (RK.     .      +RK.     .)(P.     .    ,) 
i,j+l         i,j/v    i'J  i,j-l        i,j/v    i,j-l 


IF 


irror 


-   ABSF(PN-P 

ipP- 

FAIL  =1.0 

\ 

P.     .    =    PN     - 
1,J 

CONTINUE 


L 


PRINT  FAIL 


CONTINUE 


PRINT   Output 
P.     . 


END 
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40 


42 


63 


46 


48 


43 


45 


61 


4 
13 


12 


15 


19 


PROGRAM  POTMAP 

DIMENSION  P(33,38),RK(33*38),Q(33,38) 

1,31 

,  (P(  I .J)  »J  =  1  ,12) 

1  ,31 
,  (P(  I  ,  J)  ,J  =  13,24) 

1,31 

, (P( I , J) ,  J=25,365 

1,31 

,  (Q( I , J)  »J=1»18) 

1,31 

,  (Q( I ,J) ,J=19,36) 

1,31 

,  (RK( I ,J) ,  J=l,12 ) 

1,31 

,  (RK(  I  ,J)  ,J  =  13,245 

1,31 

,  (RK(  I ,J)  ,  J=25,36) 

1H1  ) 
1,31 
,  (PC  I , J)  »J=1  ,18  ) 


1  = 
62 
1  = 
62 
1  = 
62 
1  = 
41 
1  = 
41 
1  = 
49 
1  = 
49 
1  = 
49 


DO  40 
READ 
DO  42 
READ 
DO  63 
READ 
DO  46 
READ 
DO  48 
READ 
DO  43 
READ 
DO  45 
READ 
DO  61 
READ 
PRINT  8 
FORMAT  ( 
DO   11  = 
PRINT 
PRINT 
DO   2 
PRINT 
PRINT 
DO   3 
PRINT 
PRINT 
DO   4 
PRINT 
FORMAT  ( 
LF  =  200 
ERROR  = 
DO  10  N  = 
FAIL  =  0 
DO  29  1= 
DO  29  J  = 
DO  15  K  = 
P( 10, K) = 
P( 10,14) 
PC  9,14) 
P(8,14) 
DO  19  K= 
P(7,K)  = 
PC  7,20) 
P(6,20) 
PC  5,20) 
PC  5,22)  = 
PC  5,23)  = 
PC  5,21) 
PC  5,24)  = 
PC  5,25) 
4,25) 
3,25) 
2,25) 


31 

8 

1  = 

31 

8 

1  = 

13 

8 

1  = 

13 


PC 
PC 
PC 


1,31 

9  (PC  I , J)  »J  =  19»36) 

1,31 

,  (R'<(  I  ,J),J  =  1  ,18  ) 

1,31 

, (RKC  I ,J)  ,J=19,36) 
1X,18F5.2//) 

1.0 

I9LF 
.0 

2,30 

2,35 

1,13 

PC  II9K) 

=  PC  11,15) 
=  P(9,15S 
=  P(8,15) 

15,19 

P(8,K) 
=  P(8,21S 
=  P(6,21  ) 
=  P(6,21  ) 
PC  6922) 
PC6,23  ) 
=  P(6»21  ) 
PC6924) 
=  PC6,26) 
=  PC4,26) 
=  P(3,26) 
=  PC  3*26) 


100 


1 


16 


P( 
P( 
P( 
P( 
P( 
DO 
P( 
P( 
P( 
P( 
P( 
DO 
P( 
P( 
P( 
P( 
P( 
DO 
P( 
CO 
IF 
51  IF 
520PN 
1  (  I 
2P( 
3RK 
IF 


17 


26 
27 
29 


10 

33 

30 
31 


32 
62 

49 
41 
44 


FA 
P( 
CO 
PR 

IF 
CO 
PR 
DO 
PR 
FO 
PR 
DO 
PR 
FO 
FO 
FO 
FO 
EN 


2,26) 
2,27) 
2,28) 
2,29) 
1,29) = 
18  K  = 
30, K) = 
29,18) 
30,18  ) 
28,18  ) 
27,18  ) 

16  K= 
27,  K) 
27,27) 
26,27) 
25,27) 
24,27) 

17  K  = 
24, K) 
NST  = 

(P(I  , 

(RK(  I 
=(-CON 
,J) )*P 
I,J-1) 
(  I ,J)  ) 

(ERRO 
IL  =  1 
I  ,J)  = 
NTINUE 
INT  31 

(FAIL 
NTINUE 
INT  8 

30  1  = 
INT  31 
RMAT  ( 
INT  8 

32  1  = 
INT  31 


=  P(  3 

=  P(  3 

=  P(  3 

=  P(3 

P(  1,3 

1,17 

P(29, 

=  P(29 

=  P( 

=  P( 

=  P( 

19,26 

=  P(2 

=  P( 

=  P( 

=  P( 

=  P( 

27,36 

=  P(2 

10600 

J)  )  5 

,J)  )5 

ST 

(  1-1, 
) 


,26) 
,27) 
,28  ) 

,30) 
0) 

K) 

,17) 
29,17) 
28,17) 
26,17) 

6,K) 

26,26) 

26,26) 

25,26) 

23,26) 

3,K) 

0./P( I, J) 
1,29,51 
2,29,52 

*Q(  I  »J)  +  (RK(  I 
J)  +  (RKf  T 
/  (RK(I 


+  (RK(I  +  1,J)+RK(I,J)  )*P(I+1»J)  +  (RK(I-1 
I»J+l)+RK(I,J))*P(I,J+l)+(RK(l,J-l)+RK 
1  +  1 ,J)+RK( I-1,J)+RK( I  ,J  +  i )+RK(  I ,J-1 )+4 


,J)+RK 
( I  »J)  )* 


R-ABSF(PN-P(  I  ,J)  )  )  26,27,27 
.0 
PN 


FAIL 
33,33,10 


RMAT 

RMAT( 

RMAT 

RMAT 
D 
END 


1,31 

,  (P( I, J) ,J=1,18) 

1X,18F5.0// ) 

1,31 

,  (P( I ,J) ,J=19,36 

12F6.0) 

2F6.2) 

18F4.0) 

18F4.3) 

POTMAP 


APPENDIX  IV 


FUTURE  PERFORMANCE  PREDICTION 
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FUTURE  PREDICTION  PROGRAM 
PROGRAM  PREDICT 


C 


DIMENSION   A,  B,  T,    Q,  DELP 

i 


3 


COMMON   A,  B,  TEMP,  PC;  TC,  T,  Q,  K,  DELP 


I 


) 


READ  Input  for  Z  Calc.  Constants 


X 


:> 


READ  Input  for  Q  Funct.  Tables 


ies  ] 


READ  DATA  for  Reservoir  Properti 


_L 


X 


READ  DATA  for  Past  Performance 
GAS,  DIST.,  WATER 


X 


INC  =  CPI 


r 


DO  1030    1=1,  INC,  1 


y 


READ  205  P(I) 


L 


L 


CONTINUE 


READ  Input  for  future  Performance 
GR,  DR,  WR 


Print  out  data  for  Reservoir 
Properties 
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Print  out  Past  Production  Data 

' 

B  =  PX 

1 

' 

Calculate  Z 

!' 

C  =  .00514590  (Temp) 

Bg.  =  (C)  (-) 


_  (1858064) (AK) (Time) 
D  '  (^/)(0)(Rb)Z 


I 


GRB  =  (G)(Bg.) 


FF  = 


1  -  Sgr  -  Swc 


FLOOD  =  (GRB) (FF) 


I 


Print  out  original  Conditions 


Af,  =  0 


i — ""V 


1 


DO  2    1=2,  INC 


J> 


AP.  =  P.  i  -  P. 

"   1      L-l      1 


P   =  P 

1     INC 


I 


WT  =  WATER 
DT  =  DIST 
GT  =  GAS 


W  =  (WR) (TIME) 
D  =  (DR)  (TIME)  (CONV) 
GDP  =  (GR) (TIME) 


K  =  INC 


K  =  K+l 
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L  =  K-l 


WT  =  W  +  WT 
DT  =  D  +  DT 
GT  =  GDP  +  GT 


EKAY  =  GT  +  DT 


APk  =  0.0 


Calculate  Water  Influx 


WE  =  (WIC)  (WD) 

' 

r 

POZ  =  - 

(C)G  -  EKAY) 

(GRB  +  (WT) (BW)  -  WE) 

< 

i 

Calculate  Reservoir  Pressure 

1 

1 

A  P   =  P   -  P 

A     k     1     2 

I 
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Print  Output 
Time,  Pressure,  Influx,  Gas  Produced 


Print:   All  wells 
have  flooded  out 


STOP 


>0 


Print  1  Reservoir  Depleted 


STOP 
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CALCULATION   OF  GAS   DEVIATION  FACTOR 
Subroutine   Z   CALC    (S,Z) 


(    Dii 


lmension 


s 


=    (2.) (Temp)    _ 

T 

c 

179 


PYn    =   B, 


PY2    = 
PY3   = 

PY„    = 


PY, 


(B2) (Y) 
(B3) (3Y2 
(B4) (5Y3 
(B5)(35Y^ 
(B6)(63Y" 


1) 
3Y) 

•  30Y2  +   3) 

•  70Y3  +   15Y) 


r 


X   = 


(2.)(S) 
P 


15 


14.8 


PX, 


<V 


I 


PX2   =    (B2)(X) 


PX,  = 

?\  = 

px5  = 

PX,  = 


(B3)(3X      -    1) 
(B4)(5X3    -    3X) 
(B5)(35X4    -    30X2   +   3) 
(B,)(63X5    -    70X3   +   15X) 

D 


I 


APP   =   0 

1 


DO   1 


/ 


,    M  ■-     1,6,1^ 

> 


I 


I | 


DO   1,    N  =    1,6,1 

i 


APP   =    (A        ) (PX   ) (PY   )    +  APP 
n,m  m n 


I 


Z    =   APP 


RETURN 
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CALCULATION  OF  DIMENSIONLESS  WATER  INFLUX 
Subroutine  QFUNCT  (DELT,  QT) 


QT  =  Q; 


(  Dimension  ) 

,  T  , 

(  Common  J 


TT  =  DELT 


I 


\ 


DO  3    I  =  1,381,  1  > 


DQ0  = 


Q.  -  Q.  i 
xi    i-l 

T.  -  T.  , 
l    i-l 


DQi 


T.  -  T.  . 
l    1-2 


QT  = 


Qi_2  +  [TT-Ti_2]  DQL  +  [TT-Ti_2][TT-Ti_1]  DQ2 


Continue 


Return 


End 
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WATER  INFLUX  CONDITION 
Subroutine  INFLUX  (DELTD,WD) 


C    DIMENSION  J 


f    COMMON  J 


WD  =  0.0 


mm/  DO  1    J  =  1,K,  1  \ 


TI   =   J-l 

DELT   =    (TI)  (DELTD) 

I 


Call  QFUNCT  (DELT,  QT) 


N  =  K  +  1  -  J 


I 


WI  =(DELP  )(QT) 


WD  =  WI  +  WD 


I 


CONTINUE 


RETURN 


PRESSURE  CALCULATION 
SUBROUTINE  PRESS  (Pj  POZ,  P  ) 


109 


(^DIMENSION ) 


(  COMMON  ) 


S  =  P   -  5.0 


T 


Calculate  Z 


P0ZT  =  | 


E  =  P0ZT  -  P0Z 


S=S+5.0 


S=S-5.0 


S=S+1.0 


S=S-1.0 


S=S+0.1 


S=S-0.1  [ 


S=S-.01 


>o 


S=S+.01 


no 


Ill 


800 
101 
102 


200 
1000 


201 
1010 


202 
1020 


203 


204 


205 
1030 


C 


PROGR 

DIMEN 

COMMO 

PRINT 

FORMA 

PRINT 

FORMA 

PRINT 

FORMA 

DATA 

DO  10 

DO  10 

READ 

FORMA 

CONTI 

DO  10 

READ 

FORMA 

CONTI 

DATA 

DO  10 

READ 

FORMA 

CONTI 

DATA 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

READ 

FORMA 

DATA 

READ 

FORMA 

INC  =  C 

DO  10 

READ 

FORMA 

CONTI 

DATA 


AM  PR 
SION 
N  A,B 

800 
T(  1H1 

101 
T(48H 

102 
T(  19H 
INPUT 
00  J  = 
00  1  = 
200, A 
T(  13X 
NUE 
10  1  = 
201, B 
T(  F14 
NUE 
INPUT 
20  1  = 
202, T 
T(  F12 
NUE 
INPUT 


EDICT 

A(6,6)»B(6)»T(471),Q(471)  ,DELP(999)  ,P(40) 

,TEMP,PC,TC,T,Q,K,DFLP 


RESERVOIR  COMPUTER  PROJECT FUTURE  PERFORMANCE,//) 


203 

203 

203 

203 

203 

203 

203 

203 

203 

203 

203 

203 

203 

203 

203 

203 

203 

203 

203 

203 

203 

T(  F16 

INPUT 

204, G 

T(E16 

PI 

30  1  = 

205, P 

T(F10 

NUE 

INPUT 


,A 
,C 
,R 
,D 
,T 
,S 
,0 
,A 
,S 

,u 

,A 
,T 
,S 
,S 
,P 
,P 
,T 
,B 
,C 
,W 
,G 


RESERVOIR  DATA  FOR,/) 

FOR  SUBROUTINE  ZCALC  CONSTANTS 
1,6,1 
1,6,1 
(  I  ,J) 
, F14.ll ) 

1  ,6,1 
(  I  ) 

.11  ) 

FOR  SUBROUTINE  QFUNCT  TABLE 
1 ,471,3 

(I),Q(I),T(I+l),Q(I+l),T(I+2),Q(I+2) 
.2,F12.3,F12.2,F12.3,F12.2,F12.3) 

FOR  RESERVOIR  PROPERTIES 
IWEL 
ONV 
I 

EPTH 
EMP 
GI 

H 
G 

K 

IME 

GR 

WC 

I 

C 

C 

W 

PI 

IC 

.7) 

FOR  PAST  PERFORMANCE 
AS, DIST, WATER 
.8,4X,E16.8,4X,E16.8) 

ltlNCtl 

(  I  ) 
.3) 

FOR  FUTURE  PERFORMANCE 


- 


E 


! 

i 
[ 


. 


■ 
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1  READ  700, KK 
700  FORMATf 13) 

READ  203,  GR 
READ  203,  DR 
READ  203,  WR 
C      DATA  PRINT  OUT  FOR  RESERVOIR  PROPERTIES 
PRINT  300, INC 

300  FORMATU5H  INC, CURRENT  NUMBER  OF  PRODUCTION  I  NCREMENTS=  ,  1 4  ) 
PRINT  301,AIWEL 

301  FORMAT(48H  MAX.  NUMBER  OF  PRODUCING  WELLS  AT  ANY  ONE  TIME=,F6.1) 
PRINT  302, BW 

302  FORMAT(50H  BW, FORMATION  VOLUME  FACTOR  FOR  WATER  IN  BBLS /STB= , F6 . 3 ) 
PRINT  303,CONV 

303  F0RMAT(44H  CONV,SCF  OF  VAPOR  PER  BARREL  OF  D I  ST  I LLATE= , F7. 1 ) 
PRINT  304, RI 

304  FORMAT(48H  RI, RADIUS  TO  INITIAL  EDGE  OF  RESERVOIR  IN  FEET=,F8.1) 
PRINT  305, DEPTH 

305  FORMAT(34H  DEPTH, DEPTH  OF  RESERVOIR  IN  FEET=,F8.1) 
PRINT  306, TEMP 

306  FORMAT(50H  TEMP , TEMPERATURE  OF  RESERVOIR  IN  DEGREES  RANK  I NE=  ,  F7 . 1 ) 
PRINT  307, PC 

307  FORMAT(36H  PC , PSEUDOCR I T IC AL  PRESSURE  IN  PSIA=,F6.1) 
PRINT  308, TC 

308  FORMAT(50H  TC , PSEUDOCR I T I CAL  TEMPERATURE  IN  DEGREES  RANK  I NE= , F6 . 1 ) 
PRINT  309, SGI 

309  FORMAT(40H  SGI, INITIAL  GAS  SATURATION  AS  FRACT I ON= , F5 . 3 ) 
PRINT  310,0 

310  FORMAT(24H  0, POROSITY  AS  FRACT I ON= , F5 . 3 ) 
PRINT  311, AH 

311  FORMAT(43H  AH, AVERAGE  THICKNESS  OF  RESERVOIR  IN  FEET=,F6.1) 
PRINT  312, 5G 

312  F0RMAT(44H  SG, SPECIFIC  GRAVITY  OF  GAS  RELATIVE  TO  AIR=,F7.4) 
PRINT  313, U 

313  F0RMAT(46H  U, VISCOSITY  OF  RESERVOIR  WATER  IN  CENT  I PO I SE= » F6 • 3 ) 
PRINT  314, AK 

314  FORMAT(47H  AK, WATER  PERMEABILITY  OF  RESERVOIR  IN  DARC I ES= , F7 . 4 ) 
PRINT  315, TIME 

315  FORMAT(33H  TIME, TIME  INTERVAL  USED  IN  DAYS=,F7.2) 
PRINT  316, SGR 

316  F0RMAT(41H  SGR, RESIDUAL  GAS  SATURATION  AS  FRACT I ON= , F5 . 3 ) 
PRINT  317, SWC 

317  FORMAT(42H  SWC, CONNATE  WATER  SATURATION  AS  FRACT I ON= , F5 . 3 ) 
PRINT  318, WIC 

318  FORMAT(45H  WIC, WATER  INFLUX  CONSTANT  IN  BBLS.  PER  PSI A= » F10  .  3  ,  / / ) 
C      DATA  PRINT  OUT  FOR  RESERVOIR  PRODUCTION 

PRINT  103 

103  FORMAT(55H  RESERVOIR  PRODUCTION  FOR  FUTURE  PERFORMANCE  PREDICTION) 
PRINT  104 

104  FORMATf 1H  ) 
PRINT  319, GR 

319  F0RMAT(41H  AVERAGE  DAILY  DRY  GAS  PRODUCTION  IN  SCF=,F16.1) 
PRINT  320, DR 

320  F0RMAT(44H  AVERAGE  DAILY  DISTILLATE  PRODUCTION  IN  STB=,F16.1) 
PRINT  321, WR 

321  FORMAT(39H  AVERAGE  DAILY  WATER  PRODUCTION  IN  STB= , F16 . 1 , / / ) 
C      MAIN  BODY  OF  PROGRAM 
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105 
322 
323 

400 
401 
106 
107 


402 


4 
108 

5 
6 

7 
8 


S=P( 1) 

CALL  ZCALC(S,Z) 

C=0.00514590#TEMP 

BGI=C*Z/S 

DELTD=18  58  064.0*AK*TIME/ (U*0*RI*RI ) 

GRB=G*BGI 

FF={ 1.0-SGR-SWC) /SGI 

FLOOD=GRB*FF 

PRINT  105 

FORMAT! 19H  INITIAL  CONDITIONS,/) 

PRINT  322, PI 

FORMAT(26H  INITIAL  PRESSURE  IN  PSIA=,F8.2) 

PRINT  323, G 

FORMATI29H  INITIAL  GAS  IN  PLACE  IN  SCF=,E16.8) 

PRINT  400, FLOOD 

F0RMAT(47H  MAXIMUM  FLOODABLE  RESERVOIR  VOLUME  IN  BARREL S = , E 16 . 8  ) 

PRINT  401, Z 

FORMAT(36H  INITIAL  GAS  COMPRESSIBILITY  FACTOR= , El  6 . 8 , / / ) 

PRINT  106 

F0RMAT(1H  ,/,18H  FUTURE  PREDICTION,/) 

PRINT  107 

FORMAT(15H  TIME   PRESSURE , 9X , 2HWE , 16X , 3HGPW , 8X , 7HDELTA  P) 

DELP( 1)=0.0 

DO  2  1=2, INC, 1 

DELP(  I  ) =P(  1-1 )-P( I  ) 

P1=P( INC) 

WT=WATER 

DT  =  DIST 

GT  =  GAS 

W=WR*TIME 

D=DR*TIME*CONV 

GDP=GR*TIME 

K=INC 

K  =  K  +  1 

L  =  K-1 

WT  =  W+WT 

DT=D+DT 

GT=GDP+GT 

EKAY=GT+DT 

DELP(K)=0.0 

CALL  INFLUX(DELTD,WD) 

WE=WIC*WD 

POZ=C*( G-EKAY) / (GRB+WT*RW-WE) 

CALL  PRESS(P1,P0Z,P2) 

DELP(K)=P1-P2 

PRINT  402,L,P2,WE,EKAY,DELP(K) 

FORMAT (  I5,2X,F8.2,2X,E16.8,2X,E16.8,2x,F8.2) 

P1=P2 

IF(WE-FLOOD)  8,4,4 

PRINT  108 

F0RMAT(1H  ,/,27H  ALL  WELLS  HAVE  FLOODED  OUT,/////) 

CONTINUE 


IFUK-1  ) 
CONTINUE 
GO    TO    1 
IF(P1-14.7 


6,6,7 


10,10,9 
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9 

10 

109 


GO  TO  3 
PRINT  109 
FORMAT( 1H 
CONTINUE 
GO  TO  5 
END 


/»19H  RESERVOIR  DEPLETED ,///// ) 
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SUBROUTINE  ZCA|_C(S,Z) 

DIMENSION  A(6»6)»B(6)»T(471)»Q(471)  »DELP(999  )  »PU0) 

DIMENSION  PX(6)  *PY(6) 

COMMON  A»B,TEMP,PC*TC»T ,Q,K,DELP 

Y=((2.0*TEMP/TC)-4.0)/1.9 

PY( 1) =B( 1  ) 

PY(2)=B(2 )*Y 

PY(3)=B(3)*(( 3.0*Y*Y)-1.0) 

PY(4)=B(4)*U 5.0*Y*Y)-3.0)*Y 

PY(5) =B( 5 )*(  (  (3  5.0*Y*Y ) -30 . 0 ) *Y*Y  +  3 . 0 ) 

PY(6)=B(6)*((( (63.0*Y*Y)"70.0)*Y*Y)+15.0)*Y 

X=((2.0*S/PC)-15.0)/14.8 

PX(1)=B(1) 

PX(2)=B(2)*X 

PX(3)=B(3)*( (3.0*X*X)-1.0) 

PX(4)=B(4)*( ( 5.0*X*X)-3.0)*X 

PX(5)=B(5)*((( (35.0*X*X)-30.0)*X*X)+3.0) 

PX(6)=B(6)*(  C  U  63-, 0*X*X) -70.0  )*X*X)  + 15.0) *X 

APP=0.0 

DO  1  M=l,6,l 

DO  1  N=l»6»l 

APP={ A(N,M)*PX(M)*PY(N) )+APP 

Z  =  APP 

RETURN 

END 
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SUBROUTINE  QFUNCT ( DELT » OT ) 

DIMENSION  A(6»6)>B(6)»T(471)»Q(471) ,DELP(999) »P(40) 

COMMON  A,B,TEMP,PC,TC»T,Q,K,DELP 

TT  =  DELT 

DO  3  1=1,471, 1 

IF(TT-T( I ) )  1,2,3 

DQl=(Q(I-l)-Q(I-2))/(T(I-l)-T(I-2)) 

DQ2=( (Q( I )-Q( 1-1 ) )/(T( I )-T( 1-1 ) )-DQl)/(T( I )-T( I -2) ) 

QT=Q( I-2)+(TT-T( 1-2) )*DQ1+(TT-T( 1-2) )*(TT-T( I~l ) )*DQ2 

GO  TO  4 

OT=Q( I ) 

GO  TO  4 

CONTINUE 

RETURN 

END 
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SUBROUTINE  I NFLUX ( DELTD , WD ) 

DIMENSION  A(6*6)»B(6)»T(471)»O(471)»DELP(999),P(40) 

COMMON  A»B,TEMP,PC»TC»T,Q,K,DELP 

WD=0.0 

DO  1  J=1»K,1 

TI=J-1 

DELT=TI*DELTD 

CALL  QFUNCT(DELT,QT) 

N=K+1-J 

WI=DELP(N)*QT 

WD=WI+WD 

CONTINUE 

RETURN 

END 


118 


SUBROUTINE  PRESS ( PI » POZ  ,  P2  ) 

DIMENSION  A(6  »6)  ,B(6)»T(471), 0(471)  »DELP(999)  ,P(40) 

COMMON  A,B,TEMP,PC,TC,T,Q,K,DELP 

S=Pl-5.0 

1  CALL  ZCALC(S,Z) 
POZT=S/Z 
E=POZT-POZ 
IF(ABSF(E)-5.0)  5*5*2 

2  IF(E)  4*17,3 

3  S=S-5.0 
GO  TO  1 

4  S=S+5.0 
GO  TO  1 

5  IF(ABSF(E)-1.0)  9,9*6 

6  IF(E)  8*17,7 

7  S=S-1.0 
GO  TO' 1 

8  S=S+1.0 
GO  TO  1 

9  IF(ABSF(E)-0.1 )  13,13,10 

10  IF(E)  12,17,11 

11  S=S-0.1 
GO  TO  1 

12  S=S+0.1 
GO  TO  1 

13  IF(ABSF(E)-0.02)  17,17,14 

14  IF(E)  16,17,15 

15  S=S-0.01 
GO  TO  1 

16  S=S+0.01 
GO  TO  1 

17  P2=S 
RETURN 
END 

END        PREDICT 
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DATA    INPUT    FOR    SUBROUTINE    ZCALC    CONSTANTS 


0.7071068 

1.224745 

0.7905695 

0.9354145 

0.265165 

0.293151 


2.143 

0.083 
-0.021 
•0.000 

0.004 
•0.001 

0.331 
•0.134 

0.066 
•0.027 

0.008 
•0.002 

0.105 
■0.050 

0.005 

0.010 
■0.007 

0.002 
•0.052 

0.044 
•0.019 

0.005 

0.001 
•0.002 

0.019 
■0.026 

0.019 
■0.011 

0.004 
■0.000 
■0.005 

0.008 
■0.010 

0.009 
•0.006 

0.003 


3504 

176184 

467042 

87140318 

2846283 

6595343 

23524 

03614 

880961 

174261 

8512291 

1520929 

72871 

393654 

0924798 

551336 

3181933 

6959963 

184040 

312146 

329465 

8972516 

5366676 

8326809 

703980 

383354 

262143 

535390 

2910089 

81302526 

3095900 

9178330 

894821 

559389 

0114017 

1175170 


DATA  INPUT  FOR  SUBROUTINE  QFUNCT  TABLE 


00000000 
00000010 
00000025 
00000050 
00000080 
00000200 


OOOnoooo 
00000404 
00000689 
00001020 
00001359 
00002447 


oooooooi 

00000015 
00000030 
00000060 
00000090 
00000300 


00000112 
00000520 
00000758 
00001140 
00001469 
00003202 


00000005 
00000020 
00000040 
00000070 
00000100 
00000400 


00000278 
00000606 
00000898 
00001251 
00001569 
00003893 
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00000500 
00000800 
00001100 
00001400 
00001700 
00002000 
00002300 
00002600 
00002900 
00003200 
00003500 
00003800 
00004100 
00004400 
00004700 
00005000 
00005300 
00005600 
00005900 
00006200 
00006500 
00006800 
00007100 
00007400 
00007700 
00008000 
00008300 
00008600 
00008900 
00009200 
00009500 
00009800 
00010500 
00012000 
00013500 
00015000 
00016500 
00018000 
00019500 
00021000 
00022500 
00024000 
00025500 
00027000 
00028500 
00030000 
00031500 
00033000 
00034500 
00036000 
00037500 
00039000 
00040500 
00042000 
00043500 
00045000 


^0004539 
O00O6314 
00007940 
00009461 
00010913 
00012319 
00013684 
00015013 
00016313 
00017590 
00018845 
00020080 
00021298 
00022500 
00023684 
00024855 
00026020 

00027174 
00028314 
00029443 
00030565 
00031679 
00032785 
00033883 
00034974 
00036058 
00037136 
00038207 
00039272 
00040331 
00041385 
00042433 
O0044858 
00049968 
00054976 
00059985 
00064737 
00069512 
00074226 
00078886 
00083497 
00088062 
00092589 
00097081 
00101540 
00105968 
00110367 
00114738 
00119083 
00123403 
00127699 
00131972 
00136223 
00140453 
00144664 
00148856 


00000600 
00000900 
00001200 

00001500 
00001800 
0000210Q 
00002400 
00002700 
00003000 
00003300 
00003600 
00003900 
00004200 
00004500 
00004800 
0OO05100 
00005400 
00005700 
00006000 
00006300 
00006600 
00006900 
00007200 
00007500 
00007800 
00008100 
00008400 
00008700 
00009000 
00009300 
00009600 
00009900 
00011000 
00012500 
00014000 
00015500 
00017000 
00018500 
000200^0 
OO021500 
00023000 
00024500 
00026000 
00027500 
00029000 
00030500 
00032000 
00033500 
00035000 
00036500 
00038000 
000395OH 
00041000 
00042500 
00044000 
00045500 


00005153 
0^006869 
00008457 
00009949 
00011386 
00012778 
00014131 
00015450 
00016742 
00018011 
00019259 
00020488 
00021701 
00022897 
00024076 
00025244 
00026406 
00027555 
00028691 
00029818 
00030937 
00032048 
00033151 
00034247 
00035336 
00036418 
00037494 
00038563 
00039626 
00040684 
00041735 
00042781 
00046574 
00051648 
00056625 
00061517 
00066336 
00071090 
00075785 
00080428 
00085023 
00089575 
00094090 
00098571 
00103019 
00107437 
00111827 
00116189 
00120526 
00124838 
00129126 
00133391 
00137635 
00141859 
00146064 
00150249 


00000700 
00001000 
00001300 
00001600 
00001900 
00002200 
00002500 
00002800 
00003100 
00003400 
00003700 
00004000 
00004300 
00004600 
O000490O 

O00O5200 
00005500 
00005800 
00006100 
00006400 
00006700 
00007000 
00007300 
00007600 
00007900 
00008200 
00008500 
00008800 
00009100 
00009400 
00009700 
00010000 
00011500 
00013000 
00014500 
00016000 
00017500 
0001900O 

00020500 
00022000 
00023500 
00025000 
00026500 
00028000 
00029500 
00031000 
00032500 
00034000 
00035500 
00037000 
00038500 
00040000 
00041500 
00043000 
00044500 
00046000 


00005743 
O0007411 
00008964 
00010434 
00011855 
00013233 
00014573 
00015883 
00017167 
00018429 
00019671 
00020894 
00022101 
00023291 
00024466 
00025633 
00026791 
00027935 
00029068 
00030192 
00031308 
00032417 
00033517 
00034611 
00035697 
00036777 
00037851 
00038919 
00039979 
00041034 
00042084 
00043129 
00048277 
00053317 
00058265 
00063131 
00067928 
00072661 
00077338 
00081965 
00086545 
00091084 
00095588 
00100057 
00104495 
00108904 
00113284 
00117638 
00121966 
00126270 
00130550 
00134808 
00139045 
00143262 
00147461 
00151640 
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00046500 
00048000 
00049500 
00052000 
00054000 
00057000 
00059000 
00062000 
00064000 
00067000 
00069000 
00072000 
00074000 
00077000 
00079000 
00082000 
00084000 
00087000 
00089000 
00092000 
00094000 
00097000 
00099000 
00102000 
00104000 
00107000 
00109000 
00112000 
00114000 
00117000 
00119000 
00122000 
00124000 
00127000 
00129000 
00132000 
00134000 
00137000 
00139000 
00142000 
00144000 
00147000 
00149000 
00155000 
00162500 
00170000 
00177500 
00185000 
00192500 
00200000 
00207500 
00215000 
00222500 
00230000 
00237500 
00245000 


00153029 
00157184 
00161322 
00168183 
00173639 
00181774 
00187166 
00195208 
00200542 
00208502 
00213784 
00221670 
00226904 
00234721 
00239912 
00247668 
00252819 
00260515 
00265629 
00273274 
00278353 
00285498 
00290995 
00298543 
00303560 
00311066 
00316^55 
00323517 
00328480 
00335906 
00340843 
00348230 
00353144 
00360496 
00365386 
00372704 
00377572 
00384859 
00389705 
00396959 
00401786 
00409013 
00413820 
00428196 
00446077 
00463863 
00481558 
00499167 
00516695 
00534145 
00551522 
00568830 
00586072 
00603252 
00620372 
00637437 


00047000 
00048500 
00050000 
00052500 
000550on 
00057500 
00060000 
00^625^0 
00065000 
00067500 
00070000 
00072500 
00075000 
00077500 
00080000 
00082500 
000850O0 
00087500 
00090000 
00092500 
00095000 
00097500 
00100000 
00102500 
00105000 
00107500 
001  10^00 
00112500 
00115000 
00117500 
00120000 
00122500 
00125000 
00127500 
00130000 
00132500 
00135000 
00137500 
00140000 
00142500 
00145000 
00147500 
00150000 
00157500 
00165000 
00172500 
00180000 
00187500 
00195000 
00202500 
00210000 
00217500 
00225000 
002325^0 

00240000 
002475OO 


00154416 
00158565 
00162698 
00169549 
00176357 
00183124 
00189852 
00196544 
00203201 
00209825 
00216417 
00222980 
00229514 
00236020 
00242501 
00248957 
0025  538  8 
00261795 
00268181 
00274545 
00280888 
00287211 
00293514 
00299799 
00306065 
00312314 
00318545 
00324760 
00330958 
00337142 
00343308 
00349460 
00355597 
00361720 
00367828 
00373922 
00380003 
00386070 
00392125 
00398167 
00404197 
00410214 
00416220 
00434168 
00452016 
00469771 
00487437 
00505019 
00522520 
00539945 
00557299 
00574585 
00591806 
00608965 
00626066 
00643113 


00047500 
00049000 
00051000 
00053000 
0005600^ 
00058000 
00061000 
00063000 
00066000 
00068000 
00071000 
00073000 
00076000 
00078000 
00081000 
00083000 
00086000 
00088000 
00091000 
00093000 
00096000 
00098000 
00101000 
00103000 
00106000 
00108000 
00111000 
00113000 
00116000 
00118000 
00121000 
00123000 
00126000 
00128000 
00131000 
00133000 
00136000 
00138000 
00141000 
00143000 
00146000 
00148000 
00152500 
00160000 
00167500 
00175000 
00182500 
00190000 
00197500 
00205000 
00212500 
00220000 
00227500 
00235000 
00242500 
^0250000 


00155801 
00159945 
00165444 
00170914 
00179069 
00184473 
00192533 
00197878 
00205854 
00211145 
00219046 
00224289 
00232120 
00237218 
00245086 
00250245 
00257953 
00263073 
00270729 
00275815 
00283420 
00288473 
00296030 
00301053 
00308567 
00313562 
00321032 
00326000 
00333433 
00338376 
00345770 
00350688 
00358048 
00362942 
00370267 
00375139 
00382432 
00387283 
00394543 
00399373 
00406606 
00411418 
00422214 
00440128 
00457945 
00475669 
00493307 
00510861 
00528337 
00545737 
00563068 
00580332 
00597532 
00614672 
00631755 
00648781 
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00255000 
00270000 
00285000 
00300000 
00315000 
00330000 
00345000 
00360000 
00375000 
00390000 
00405000 
00420000 
00435000 
00450000 
00465000 
00480000 
00495000 
00520000 
00550000 
00580000 
00610000 
00640000 
00670000 
00700000 
00730000 
00760000 
00790000 
00820000 
00850000 
00880000 
00910000 
00940000 
00970000 
01000000 
01750000 
03000000 
05000000 
07500000 
10000000 


00660093 
00693877 
00727449 
00760833 
00794042 
00827088 
00859974 
00892712 
00925309 
00957773 
00990108 
01022318 
01054409 
01086390 
01118257 
01150012 
01181666 
01234203 
01296893 
01359225 
01421224 
01482912 
01544305 
01605418 
01666265 
01726859 
01787212 
01847336 
01907243 
01966943 
02026438 
02085744 
02144878 
02203861 
03633368 
05899508 
09363099 
13531457 
17586284 


00260000 
00275000 
00290000 
00305000 
00320000 
00335000 
003500^0 
003650H0 
00380000 
00395000 
00410000 
00425000 
00440000 
00455000 
00470000 
00485000 
00500000 
00530000 
00560000 
00590000 
00620000 
00650000 
00680000 
00710000 
00740000 
00770000 
00800000 
00830000 
00860000 
00890000 
00920000 
00950000 
00980000 
01250000 
020000O0 
03500000 
06000000 
08000000 
12500000 


00671379 
00705090 
00738598 
00771922 
00805075 
00838067 
00870903 
00903594 
00936144 
00968566 
01000858 
01033028 
01065082 
01097024 
01128854 
01160574 
01192198 
01255141 
01317709 
01379927 
01441820 
01503408 
01564706 
01625729 
01686490 
01747002 
01807278 
01867329 
01927166 
01986796 
02046227 
02105473 
02164555 
02688967 
04095800 
06780247 
11047299 
14350121 
21560732 


00265000 
00280000 
00295000 
00310000 
00325000 
00340000 
00355000 

00370000 
00385000 
00400000 
00415000 
00430000 
00445000 
00460000 
00475000 
00490000 
00510000 
00540000 
00570000 
00600000 
00630000 
00660000 
00690000 
00720000 
00750000 
00780000 
00810000 
00840000 
00870000 
00900000 
00930000 
00960000 
00990000 
01500000 
02500000 
04000000 
07000000 
09000000 
15000000 


00682640 
00716280 
00749725 
00782992 
00816090 
00849028 
00881816 
00914459 
00946966 
00979344 
01011595 
01043724 
01075743 
01107646 
01139439 
01171125 
01213222 
01276037 
01338486 
01400593 
01462383 
01523872 
01585077 
01646011 
01706688 
01767120 
01827319 
01887298 
01947065 
02006628 
02065996 
02125184 
02184216 
03164780 
05005726 
07650096 
12708358 
15975389 
25380000 


C 


DATA  INPUT  FOR  RESERVOIR  PROPERTIES 


2.0 

768.0 

2190. 

11450. 

690. 

.80 

.28 

64.4 

.7975 

.30 

.050 

180. 

.30 
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0.  15 

8838. 

662.26 

387.48 

1.03 

7.0 

7.52757 

606898622000. 

2.07907540E+10 
8838. 
8764. 
8657. 
8509. 
8364. 
8219. 
8100. 

9 
25000000. 
1250. 
30. 


1.54221800E+06 


8.67600000F+03 


RESERVOIR  COMPUTER  PROJECT FUTURE  PERFORMANCE 


RESERVOIR  DATA  FOR 


124 


INC.C 

MAX. 

BW,FO 

CONV, 

RI  ,RA 

DEPTH 

TEMP, 

PC, PS 

TCPS 

SGI  ,1 

0,POR 

AH,AV 

SG,SP 

U,VIS 

AK,WA 

TIME, 

SGR,R 

SWCC 

WICW 


URREN 

NUMBE 

RMATI 

SCF  0 

DIUS 

,DEPT 

TEMPE 

EUDOC 

EUDOC 

NITIA 

OSITY 

ERAGE 

ECIFI 

COSIT 

TER  P 

TIME 

ESIDU 

ONNAT 

ATER 


T  NUMBER 
R  OF  PRO 
ON  VOLUM 
F  VAPOR 
TO  INITI 
H  OF  RES 
RATURE  0 
RITICAL 
RI TICAL 
L  GAS  SA 
AS  FRAC 
THICKNE 
C  GRAVIT 
Y  OF  RES 
FRMEARIL 
INTERVAL 
AL  GAS  S 
E  WATER 
INFLUX  C 


OF  P 
DUCIN 
E  FAC 
PER  B 
AL  ED 
ERVOI 
F  RES 
PRESS 
TEMPE 
TURAT 
TI0N  = 
SS  OF 
Y  OF 
ERVOI 
ITY  0 

USED 
ATURA 
SATUR 
ONSTA 


RODUCTIO 
G  WELLS 
TOR  FOR 
ARREL  OF 
GE  OF  RE 
R  IN  FEE 


ERVOIR 
URfe  IN 
RATURE 
ION  AS 

.280 

RESERVO 
GAS  RELA 
R  WATER 
F  RESFRV 

TN  DAYS 
TION  AS 
ATION  AS 
NT  IN  BB 


N  INCREMENTS=    7 

AT  ANY  ONE  TIME=    2.0 

WATER  IN  BBLS/STB=  1.030 

DISTILLATE=   768.0 
SERVOIR  IN  FEET=   2190.0 
T=  11450.0 

N  DEGREES  RANKINE=   690.0 
SIA=  662.3 

N  DEGREES  RANKINE=  387.5 
RACTI0N=  .800 

IR  IN  FEET=   64.4 
TIVE  TO  AIR=   .7975 
IN  CENTIP0ISE=   .300 
OIR  IN  DARCTES=   .0500 
=  180.00 
FRACTION=  .300 

FRACTI0N=  .150 
LS.  PER  PSIA=      7.528 


RESERVOIR  PRODUCTION  FOR  FUTURE  PERFORMANCE  PREDICTION 

AVERAGE  DAILY  DRY  GAS  PRODUCTION  IN  SCF=  25000000.0 

AVERAGE  DAILY  DISTILLATE  PRODUCTION  IN  STB=  1250.0 

AVERAGE  DAILY  WATER  PRODUCTION  IN  STB=  30.0 


INITIAL  CONDITIONS 


INITIAL  PRESSURE  IN  PSIA=  8838.00 
INITIAL  GAS  IN  PLACE  IN  SCF=   6.06898622E  11 
MAXIMUM  FLOODABLE  RESERVOIR  VOLUME  IN  BARRELS= 
INITIAL  GAS  COMPRESSIBILITY  FACTOR=   1.32981926E 


2.22914288E 
00 


08 


FUTURE  PREDICTION 


TIME   PRESSURE 


WE 


GPW 


DELTA  P 


7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 


8009.66 
7870.52 
7734.90 
7602.48 
7473.08 
7346.57 
7222.84 
7101.81 
6983.39 
6867.53 
6754.  18 


3.019976 
3.927644 
4.987118 
172960 
477339 
894323 
041889 
1.204660 
1.377350 
1.559595 
1.751055 


62E 

87E 
13E 
89E 
89E 
HE 
07E 
79E 
30E 
84E 
77E 


05 
05 
05 
05 
05 
05 
06 
06 
06 
06 
06 


2.546509 
3.013789 
3.481069 
3.948349 
4.415629 
4.882909 
5.350189 
5.817469 
6.284749 
6.752029 
7.219309 


62F 
62E 
62E 
62E 
62E 
62E 
62E 
62E 
62E 
62F 
62E 


10 
10 
10 
10 
10 
10 
10 
10 
10 
10 
10 


90.34 
139.14 
135.62 
132.42 
129.40 
126.51 


123 
121 
118 


73 
03 
42 


115.86 
113.35 
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18 

6643.27 

1.95140676E 

06 

7.68658962E 

10 

110.91 

19 

6534.76 

2.16034859E 

06 

8.  15386962E 

10 

108.51 

20 

6428.60 

2.37757163E 

06 

8.62114962E 

10 

106.16 

21 

6324.75 

2.60283387E 

06 

9.08842962E 

10 

103.85 

22 

6223.  15 

2.83585730E 

06 

9.55570962E 

10 

101.60 

23 

6123.75 

3.07640675E 

06 

1.00229896E 

11 

99.40 

24 

6026.50 

3.32422979E 

06 

1.04902696E 

11 

97.25 

25 

5931.35 

3.57910580E 

06 

1.09575496E 

11 

95.15 

26 

5838.24 

3.84080861E 

06 

1.14248296E 

11 

93.11 

27 

5747.13 

4.10910827E 

06 

1.18921096E 

11 

91.11 

28 

5657.94 

4.38383777E 

06 

1.23593896E 

11 

89.19 

29 

5570.63 

4.66479429E 

06 

1.28266696E 

11 

87.31 

30 

5485.  14 

4.95179404E 

06 

1.32939496E 

11 

85.49 

31 

5401.40 

5.24466308E 

06 

1.37612296E 

11 

83.74 

32 

5319.37 

5.54323484E 

06 

1.42285096E 

11 

82.03 

33 

5238.99 

5.84734876E 

06 

1.46957896E 

11 

80.38 

34 

5160.  19 

6.15685365E 

06 

1.51630696E 

11 

78.80 

35 

5082.93 

6.47160388E 

06 

1.56303496E 

11 

77.26 

36 

5007.  15 

6.79146218E 

06 

1.60976296E 

11 

75.78 

37 

4932.79 

7.11629632E 

06 

1.65649096E 

11 

74.36 

38 

4859.81 

7.44598239E 

06 

1.70321896E 

11 

72.98 

39 

4788.  15 

7.78039908E 

06 

1.74994696E 

11 

71.66 

40 

4717.76 

8.  11943332E 

06 

1.79667496E 

11 

70.39 

41 

4648.59 

8.46297811E 

06 

1.84340296E 

11 

69.17 

42 

4580.60 

8.81092942E 

06 

1.89013096E 

11 

67.99 

43 

4513.74 

9.16318972E 

06 

1.93685896E 

11 

66.86 

44 

4447.97 

9.51966598E 

06 

1.98358696E 

11 

65.77 

45 

4383.24 

9.88026859E 

06 

2.03031496E 

11 

64.73 

46 

4319.51 

1.02449144E 

07 

2.07704296E 

11 

63.73 

47 

4256.74 

1.06135227E 

07 

2.12377096E 

11 

62.77 

48 

4194.90 

1.09860180E 

07 

2.17049896E 

11 

61.84 

49 

4133.94 

1.13623269E 

07 

2.21722696E 

11 

60.96 

50 

4073.83 

1.17423812E 

07 

2.26395496E 

11 

60.11 

51 

4014.54 

1.21261164E 

07 

2.31068296E 

11 

59.29 

52 

3956.03 

1.25134699E 

07 

2.35741096E 

11 

58.51 

53 

3898.27 

1.29043838E 

07 

2.40413896E 

11 

57.76 

54 

3841.23 

1.32988022E 

07 

2.45086696E 

11 

57.04 

55 

3784.88 

1.36966727E 

07 

2.49759496E 

11 

56.35 

56 

3729.  18 

1.40979449E 

07 

2.54432296E 

11 

55.70 

57 

3674.  12 

1.45025730E 

07 

2.59105096E 

11 

55.06 

58 

3619.67 

1.49105112E 

07 

2.63777896E 

11 

54.45 

59 

3565.79 

1.53217166E 

07 

2.68450696E 

11 

53.88 

60 

3512.46 

1.57361517E 

07 

2.73123496E 

11 

53.33 

61 

3459.66 

1.61537783E 

07 

2.77796296E 

11 

52.80 

62 

3407.38 

1.65745573E 

07 

2.82469096E 

11 

52.28 

63 

3355.57 

1.69984484E 

07 

2.87141896E 

11 

51.81 

64 

3304.23 

1.74254188E 

07 

2.91814696E 

11 

51.34 

65 

3253.32 

1.78554343E 

07 

2.96487496E 

11 

50.91 

66 

3202.82 

1.82884669E 

07 

3.01160296E 

11 

50.50 

67 

3152.74 

1.87244912E 

07 

3.05833096E 

11 

50.08 

68 

3103.03 

1.91634781E 

07 

3.  10505896E 

11 

49.71 

69 

3053.68 

1.96054041E 

07 

3.15178696E 

11 

49.35 

70 

3004.67 

2.00502457E 

07 

3.19851496E 

11 

49.01 

71 

2955.97 

2.04979827E 

07 

3.24524296E 

11 

48.70 

72 

2907.60 

2.09485976E 

07 

3.29197096E 

11 

48.37 

73 

2859.50 

2.14020682E 

07 

3.33869896E 

11 

48.10 
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74 

2811. 

.68 

2.18583818E 

07 

3.38542696E 

11 

75 

2764, 

,11 

2.23175214E 

07 

3.43215496E 

11 

76 

2716. 

,78 

2.27794740E 

07 

3.47888296E 

11 

77 

2669, 

,68 

2.32442265E 

07 

3.52561096E 

11 

78 

2622, 

,78 

2.37117672E 

07 

3.57233896E 

11 

79 

2576, 

,07 

2.41820867E 

07 

3.61906696E 

11 

80 

2529, 

,54 

2.46551770E 

07 

3.66579496E 

11 

81 

2483, 

,18 

2.51310297E 

07 

3.71252296E 

11 

82 

2436, 

,96 

2.56096366E 

07 

3.75925096E 

11 

83 

2390, 

,86 

2.60909930E 

07 

3.80597896E 

11 

84 

2344. 

,88 

2.65750960E 

07 

3.85270696E 

11 

85 

2299. 

,05 

2.70619420E 

07 

3.89943496E 

11 

86 

2253. 

,25 

2.75515216E 

07 

3.94616296E 

11 

87 

2207. 

,58 

2.80438455E 

07 

3.99289096E 

11 

88 

2161, 

,95 

2.85389026E 

07 

4.03961896E 

11 

89 

2116. 

,35 

2.90366994E 

07 

4.08634696E 

11 

90 

2070, 

,82 

2.95372392E 

07 

4.13307496E 

11 

91 

2025, 

,29 

3.00405193E 

07 

4.17980296E 

11 

92 

1979. 

,77 

3.05465482E 

07 

4.22653096E 

11 

93 

1934. 

,24 

3.10553302E 

07 

4.27325896E 

11 

94 

1888. 

,69 

3.15668727E 

07 

4.31998696E 

11 

95 

1843. 

,09 

3.20811827E 

07 

4.36671496E 

11 

96 

1797, 

,46 

3.25982721E 

07 

4.41344296E 

11 

97 

1751. 

,76 

3.31181457E 

07 

4.46017096E 

11 

98 

1705, 

,96 

3.36408157E 

07 

4.50689896E 

11 

99 

1660, 

,09 

3.41662966E 

07 

4.55362696E 

11 

100 

1614. 

,10 

3.46945965E 

07 

4.60035496E 

11 

101 

1568. 

,00 

3.52257331E 

07 

4.64708296E 

11 

102 

1521. 

,77 

3.57597186E 

07 

4.69381096E 

11 

103 

1475. 

,36 

3.62965691E 

07 

4.74053896E 

11 

104 

1428. 

,76 

3.68363065E 

07 

4.78726696E 

11 

105 

1382. 

01 

3.73789519E 

07 

4.83399496E 

11 

106 

1335, 

01 

3.79245192E 

07 

4.88072296E 

11 

107 

1287, 

,81 

3.84730397E 

07 

4.92745096E 

11 

108 

1240, 

,39 

3.90245323E 

07 

4.97417896E 

11 

109 

1192, 

,68 

3.95790186E 

07 

5.02090696E 

11 

110 

1144. 

,68 

4.01365305E 

07 

5.06763496E 

11 

111 

1096. 

,38 

4.06970957E 

07 

5.11436296E 

11 

112 

11047. 

,78 

4.12607427E 

07 

5.16109096E 

11 

113 

998. 

79 

4.18275004E 

07 

5.20781896E 

11 

114 

949. 

48 

4.23974101E 

07 

5.25454696E 

11 

115 

899. 

74 

4.29704978E 

07 

5.30127496E 

11 

116 

849. 

62 

4.35468084E 

07 

5.34800296E 

11 

117 

799. 

02 

4.41263736E 

07 

5.39473096E 

11 

118 

747. 

99 

4.47092427E 

07 

5.44145896E 

11 

119 

696. 

45 

4.52954528E 

07 

5.48818696E 

11 

120 

644. 

,35 

4.58850553E 

07 

5.53491496E 

11 

121 

591. 

74 

4.64781051E 

07 

5.58164296E 

11 

122 

538. 

,51 

4.70746455E 

07 

5.62837096E 

11 

123 

484. 

,66 

4.76747371E 

07 

5.67509896E 

11 

124 

430. 

,14 

4.82784344E 

07 

5.72182696E 

11 

125 

374. 

,92 

4.88857997E 

07 

5.76855496E 

11 

126 

318. 

,93 

4.94968966E 

07 

5.81528296E 

11 

127 

262. 

,20 

5.01117984E 

07 

5.86201096E 

11 

128 

204. 

,60 

5.07305695E 

07 

5.90873896E 

11 

129 

146. 

,  12 

5.13532948E 

07 

5.95546696E 

11 

47.82 

47.57 
47.33 
47.10 
46.90 
46.71 
46.53 
46.36 
46.22 
46.10 
45.98 
45.83 
45.80 
45.67 
45.63 
45.60 
45.53 
45.53 
45.52 
45.53 
45.55 
45.60 
45.63 
45.70 
45.80 
45.87 
45.99 
46.10 
46.23 
46.41 
46.60 
46.75 
47.00 
47.20 
47.42 
47.71 
48.00 
48.30 
48.60 
48.99 
49.31 
49.74 
50.12 
50.60 
51.03 
51.54 
52.10 
52.61 
53.23 
53.85 
54.52 
55.22 
55.99 
56.73 
57.60 
58.48 


130 

86.70 

131 

26.28 

132 

-35.20 

5.19800526E  07 
5.26109289E  07 
5.32460162E  07 


6.00219496E  11  59. 42 
6.04892296E  11  60.42 
6.09565096E  11      61.48 
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VITA 


Ossian  R.  Butterfield  was  born  at  Leominster,  Massa- 
chusetts, on  May  12,  1927,  the  son  of  Ossian  R.  and  Helen  P. 
Butterfield.   He  attended  public  schools  in  Leominster, 
graduating  from  Leominster  High  School  in  1944.   He  entered 
the  U.  S.  Navy  in  June,  1944,  was  enrolled  in  the  Navy  V-12 
College  Training  Program  and  attended  Williams  College  in 
W i  1  1  iams town ,  Massachusetts.   In  November,  1945,  he  was 
transferred  to  the  U.S.  Naval  Reserve  Officers  Training 
Course  Unit  at  Brown  University.    In  1947,  he  was  graduated 
from  Brown  with  the  degree  of  Bachelor  of  Science  in  Engi- 
neering.  Upon  graduation,  he  was  commissioned  an  Ensign  in 
the  U.S.  Navy  Civil  Engineering  Corps.   Duty  assignments  have 
included  engineering  duties  at  the  Naval  Proving  Grounds, 
Dahlgren,  Virginia;  the  Naval  Air  Station,  Dallas,  Texas;  the 
Naval  Stations  in  Argentia,  Newfoundland,  and  in  Bermuda,  The 
West  Indies;  the  Naval  Air  Station,  Patuxent  River,  Maryland; 
the  Naval  Radio  Station,  Kamiseya,  Japan;  and  the  Naval  Mine 
Warfare  School,  Yorktown,  Virginia.    In  1955  and  1956  he 
attended  Rensselaer  Polytechnic  Institute  in  Troy,  New  York, 
where  he  received  the  Bachelor  of  Civil  Engineering  degree  in 
June  1956.   He  has  served  as  Executive  Officer  and  as  Command- 
ing Officer  of  U.S.  Naval  Mobile  Construction  Battalion  Six 
deployed  to  Rota,  Spain,  and  other  bases  in  the  Atlantic  Ocean 


Area.   His  most  recent  duty  assignment  was  in  the  Bureau  of 
Yards  and  Docks,  Navy  Department,  Washington,  D.C.,  where  he 
was  Director  of  the  Construction  Division  responsible  for  the 
management  of  projects  engineered  and  constructed  by  the  U.S. 
Navy  for  other  governmental  agencies. 

He  was  selected  for  postgraduate  education  in  petroleum 
engineering  and  entered  the  Graduate  School  of  The  University 
of  Texas  in  September,  1963. 

Butterfield  was  married  to  Ann  Churchill  of  King  George, 
Virginia,  on  April  11,  1950.   They  have  two  children,  Marcia, 
age  13,  and  Brian,  age  12. 

Permanent  address:   6600  Willamette  Drive 

Aus  t in ,  Texas 


This  thesis  was  typed  by  Virginia  F.  Dailey, 


DUDLEY  KNOX  LIBRARY 


3  2768  00032593  0 


